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ABSTRACT 


A  program  of  investigation  was  begun  in  April  1964,  to  determine 
nondestructive  testing  methods  and  techniques  for  the  evaluation  and 
characterization  of  graphite  materials  under  contract  No.  AF33(615)-1601. 

During  the  first  year,  April  1964  to  April  1965,  the  properties  and  behavior 
characteristics  which  are  necessary  for  characterizing  graphite  as  an 
ablative  material  for  aerospace  applications  were  identified,  and  applicable 
NDT  methods  and  techniques  were  described  and  theoretically  justified. 
Moreover,  the  philosophy  of  Avco's  step-wise  approach  to  nondestructive 
materials  evaluation  was  identified  in  relation  to  graphite,  and  an  over-all 
program  plan  was  established.  Correlations  between  NDT  measurements 
and  properties  of  interest  were  developed  and  reported  for  certain  grades 
of  graphite.  During  the  second  year,  April  1965  to  April  1966,  these 
correlations  were  extended  to  other  grades  and,  additionally,  have  been 
studied  in  combination  in  an  effort  to  obtain  qualitative  response  to 
behavior  characteristics  such  as  erosion  rate  and  thermal  shock  resistance. 

To  this  end,  nondestructive  test  methods  for  evaluating  the  performance 
of  graphite  under  complex  thermal  and  mechanical  loading  conditions  were 
investigated.  The  observed  response  of  a  molded  graphite  is  compared 
to  theoretical  predictions,  which  are  based  on  the  constitutive  equations 
for  a  transversely  isotropic  media.  Also  discussed  are  the  limitations 
of  the  theory  of  anisotropic  elasticity  as  applied  to  the  graphite.  Theoretical 
implications  of  thermal  stress  wave  phenomena  are  discussed,  and  an 
attempt  is  made  to  gain  insight  into  the  failure  mechanisms  of  such  materials. 

It  is  concluded  that  a  temperature  dependent  Mohr  Coulomb  yield  envelope 
may  be  adequate  to  predict  certain  aspects  of  failure  in  graphite.  In  addition, 
realistic  thermal  shock  parameters  might  be  formulated  by  using  detailed 
experimental  investigations  of  wave  phenomena  as  well  as  asymptotic  solutions 
to  the  coupled  anisotropic  thermoelastic  field  equations. 

In  connection  with  nondestructive  thermal  properties  measurements  and 
related  thermal  shock  studies,  a  novel  technique  has  been  developed  to 
nondestructively  determine  the  near- room-temperature  thermal  conductivity 
of  graphite.  Theory  of  the  technique  is  presented  and  is  supported  by  limited 
experimentation  on  five  good  conductors,  with  emphasis  on  grain-oriented 
graphites,  to  demonstrate  its  feasibility.  The  technique  involves  simultaneously 
heating  the  surface  of  a  "semi-infinite"  solid  with  radiant  energy  and  obaerving 
the  temperature  history  of  this  surface  using  an  infrared  radiometer  as  the 
radiance  serior.  Analysis  of  this  history,  together  with  knowledge  of  the  solid's 
density  and  heat  capacity  and  comparison  with  a  standard  history,  yields 
the  conductivity.  Results  indicate  the  accuracy  of  the  technique  to  be 
within  17  percent  of  the  standard  using  present  equipment. 
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SECTION  I. 


INTRODUCTION 


Material  design  parameters  and  related  properties  requirements  are 
among  the  more  important  considerations  in  the  design  of  aerospace  hardware. 
Inconsistency  in  the  behavior  of  graphite  materials  in  ablative  applications  can 
be  traced  to  a  significant  degree  to  the  lack  of  adequate  knowledge  of  material 
propex.^es  in  relation  to  performance  data.  The  inherent  variability  within  and 
between  billets  of  various  aerospace  grades  of  graphite  is  evident  from  the 
widespread  properties  data  which  have  been  reported.  Such  evidence  suggests 
that  although  ineffective  design  models  for  ablative  applications  are  a  contributing 
factor,  the  material  variability  is  the  predominant  influence  affecting  graphite's 
inconsistent  behavior.  Moreover,  because  of  the  wide  range  in  properties  that 
can  exist  in  graphite  statistical  analysis  of  random  destructive  data  for  a  given 
grade  is  not  sufficient  to  characterize  a  given  billet  or  hardware  component  for 
ablative  environments.  As  a  result,  interest  has  turned  to  nondestructive  testing 
(NDT)  techniques  for  characterizing  graphite  for  these  applications. 

The  NDT  evaluation  program  which  has  been  conducted  on  graphite  at  Avco/RAD 
was  begun  in  April  1964.  The  results  of  the  first  year's  work  were  published  in 
the  final  report  for  that  period.  ^  The  over-all  objective  of  the  work  has  been  to 
develop  meaningful  nondestructive  tests  for  the  evaluation  of  graphite  in  thicknesses 
up  to  5  inches.  A  major  part  of  this  activity,  which  has  carried  over  into  the 
second  year,  has  been  concerned  with  the  development  of  correlations  between 
destructive  and  nondestructive  tests  for  significant  design  properties.  The 
success  of  Avco/RAD's  approach  in  this  effort  has  been  described  in  detail  in 
reference  1.  Briefly,  this  approach  depends  upon  the  following: 

l)Accurate  definition  of  the  service  significant  properties,  including  a 
knowledge  of  the  r  roblem  areas  related  to  performance  and  the  failure 
modes  and  mechanisms  which  predominate. 

Z)  A  knowledge  of  the  physical  composition  o+  the  material,  and  the 
variables  influencing  the  properties  of  interest. 

3)  Selection  of  the  NDT  techniques  which  theoretically  should  exhibit  an 
interaction  with  these  material  variables  and,  hence,  respond  to  the 
associated  properties. 

4)  Experimental  verification  of  the  selected  NDT  techniques  by  correlating 
destructively  determined  properties  with  NDT  measurements. 

This  report  discusses  the  work  performed  on  this  contract  during  the  second 
year,  April  1965  to  April  1966.  During  this  period,  correlations  between  NDT 
measurements  and  properties  of  interest  have  been  extended  to  include 


all  the  grades  identified  to  be  of  interest  to  the  Air  Force.  Additionally,  combined 
use  of  various  NDT  techniques  and  related  correlations  have  been  investigated  in 
an  effort  to  obtain  qualitative  response  to  behavior  characteristics  such  as 
erosion  rate  and  thermal  shock  resistance.  Specifically,  NDT  methods  and 
techniques  for  evaluating  the  performance  of  graphite  under  complex  thermal 
and  mechanical  loading  conditions  were  investigated,  and  a  novql  technique 
for  nondestructively  determining  the  near-room-temperature  thermal  conductivity 
of  graphite  has  been  developed.  Additional  emphasis  during  the  past  year  has  been 
devoted  to  the  selection  of  devices  which  can  readily  be  applied  to  practical 
inspection  problems  concerned  with  the  characterization  of  graphite  for 
ablative  applications. 
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SECTION  II. 

LITERATURE  AND  INDUSTRIAL  SURVEY 


The  literature  and  consultation  survey,  which  was  begun  in  April  1964,  has  been 
carried  out  as  a  continuing  part  of  the  over-all  program.  Initially,  this  survey  was 
aimed  at  obtaining  information  concerning  destructive  and  nondestructive  testing  of 
graphite,  definition  and  ranking  of  the  service-significant  properties  and  behavior 
characteristics,  and  the  effects  of  material  processing  variables  on  these  service 
requirements.  During  the  past  year,  information  was  also  sought  as  to  the 
performance  of  graphite  components  during  environmental  testing,  including  failure 
modes  and  mechanisms  in  an  effort  to  determine  the  extent  to  which  such  failures 
could  be  identified  with  materials  variability.  To  date,  however,  only  limited 
information  has  been  revealed  which  attempts  to  link  graphite  material  properties 
variations  to  environmental  performance,  and  this  information  is  mostly  concerned 
with  erosion  rate  for  subscale  rocket  nozzle  inserts  as  related  to  porosity. 


SECTION  III 


COMPILATION  OF  EXPERIMENTAL  AND  THEORETICAL  DATA 


The  experimental  and  theoretical  compilation  has  also  continued  throughout 
the  second  year.  A  detailed  discussion  of  this  work  for  the  first  year,  including 
tables  relating  property  and  behavior  characteristics  to  influencing  materials 
variables,  and  the  interaction  of  the  selected  NDT  energy  with  these  variables, 
was  also  presented  in  Reference  1.  Also  included  in  that  discussion  was  a 
theoretical  justification  for  applying  each  selected  NDT  technique,  namely, 
ultrasonic  velocity  and  its  relation  to  modulus  and  density,  radiography  for 
flaw  detection,  and  radiometry  for  measurement  of  local  bulk  density.  Theory 
and  justification  for  use  of  transient  infrared  techniques  for  determining  thermal 
properties  were  also  presented  briefly,  and  will  be  covered  in  detail  in  a  later 
section  of  this  report.  Each  of  these  NDT  techniques  is  applicable  to  practical 
inspection  and  evaluation  of  graphite  billets  and  hardware  components. 

From  the  more  significant  design  properties  for  ablative  applications, 
modulus  of  elasticity,  tensile  strength  and  strain-to-failure,  thermal  conductivity 
and  thermal  expansion  were  selected  for  study.  Developing  the  NDT  capability 
for  determining  each  of  these  properties  at  room  temperature,  or  near-room 
temperature,  is  justification  for  the  correlation  studies  which  have  been  performed 
on  each  property.  However,  in  addition  to  this  single  importance  of  each  of  these 
properties  is  their  combined  significance  in  relation  to  behavior  characteristics 
such  as  erosion  rate  and  thermal  shock  resistance. 


The  importance  of  erosion  rate  and  thermal  shock  resistance  in  graphite 
rocket  nozzle  applications,  has  been  noted  by  other  investigators.  '  ’  ' 

These  workers  report  on  the  effects  of  porosity  on  erosion  rate  for  subscale 
nozzle  inserts  during  solid  propellant  rocket  motor  firings.  C.  E.  Waylett 
et  al  (Reference  2),  for  example,  show  significant  data  for  an  erosion  rate 
parameter  versus  bulk  density  for  nozzle  inserts  of  several  grades  of  graphite. 
This  is  particularly  significant  from  the  viewpoint  of  NDT  since  local  variations 
in  bulk  density  of  discrete  volumes  in  billets  and  nozzles  can  be  determined 
nondestructively  by  radiometric  gauging.  Berard  and  Swope  (Reference  3) 
describe  the  effects  of  size  and  distribution  of  open  porosity  on  erosion  rate 
for  some  fifteen  grades  of  graphite.  Moreover,  they  give  a  detailed  description, 
as  do  other  sources  in  the  literature,  ^  of  how  to  measure  the  size  of  open 
porosity  by  mercury  intrusion  techniques.  Again,  from  the  viewpoint  of  NDT, 
measurement  of  open  porosity  coupled  with  local  measurements  of  bulk  density, 
from  which  distribution  of  total  volume  porosity  can  be  assessed,  is  significant. 


The  occurrence  of  failures  caused  by  thermal  shock  and  related  properties 
variability  are  not  reported  by  either  of  the  investigating  teams  referenced 
above,  although  the  ability  to  withstand  thermal  shock  is  recognized  by  Waylett, 
et  al,  as  being  a  prime  factor  for  determining  tne  suitability  of  graphite  for 
nozzle  applications.  It  should  be  mentioned  that,  since  neither  of  the  programs  referenced 
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above  was  designed  to  evaluate  thermal  shock  resistance,  it  may  be  that  the 
environmental  conditions  used  were  not  sufficient  in  degree  and/or  duration 
to  induce  this  kind  of  failure. 


In  any  case,  data  on  ablative  performance  of  graphite  materials  versus 
thermal  shock  and  related  properties  has  not  been  noted  in  the  literature  to 
date.  However,  references  to  a  qualitative  parameter  for  describing  the 
resistance  of  high-temperature  refractory  materials  to  fracture  by  thermal 
shock  have  been  noted  in  the  literature  ^).  This  parameter  has  been 
expressed  as: 


F 


kS 
a  E 


where:  k 

S 
a 
E 


thermal  conductivity 
tensile  strength 

coefficient  of  thermal  expansion 
modulus  of  elasticity 


(1) 


Nondestructive  determinations  of  tensile  strength  and  modulus  have  been  presented 
previously,  and  will  be  discussed  in  some  detail  in  Section  IV  of  this  report. 

The  coefficient  of  thermal  expansion  and  thermal  conductivity,  however,  are  not 
easily  handled  insofar  as  nondestructive  measurement  is  concerned,  although 
the  possibility  of  obtaining  a  nondestructive  assessment  of  thermal  expansion 
is  suggested  from  the  following  relationship  between  thermal  expansion  and 
other  properties  in  crystalline  isotropic  metals: 


K  a 

-  =  Q(a  constant) 

pcv 


where: 


K 


a 


P 


bulk  modulus 

coefficient  of  thermal  expansion 
density 

heat  capacity  at  constant  volume 


Substituting  for  a  from  equation  (2)  into  equation  (1),  we  obtain: 


(2) 


F 


k  KS 

pCv  EQ 


(3) 
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D.  Preliminary  Investigation  of  the  Thermo- Mechanical  Behavior  of  Graphite 

(7  8) 

Several  investigators  ’  have  suggested  that,  because  of  the  complex 
combinations  of  failure  mechanisms  that  may  occur,  prediction  of  the 
behavior  of  potential  rocket  nozzle  materials  based  on  physical  property 
data  or  simple  laboratory  tests  is  difficult,  if  not  impossible. 

For  graphite,  this  is  in  part  due  to  ignorance  of  the  true  failure  mechanisms 
of  ouch  materials  and  also  partly  attributable  to  a  general  over -confidence 
in  simple  mathematical  idealization  sf  graphitic  mechanical  and  thermal 
behavior.  In  this  discussion,  the  widely  applied  assumption  of  plane  isotropy 
is  first  presented,  and  an  effort  is  then  made  to  delineate  some  of  the 
controlling  influences  of  actual  response  phenomena. 

Typical  microstructure  of  molded  ATJ  graphite  is  illustrated  by  the 
photomicrographs,  Figures  8  and  9.  The  non-uniform  density  as  well  as 
the  inherent  porosity  is  clearly  revealed  by  these  magnified  views  of  the 
material.  Due  to  observed  anisotropy  and  the  nature  of  the  molding  process, 
a  widely  suggested  form  of  constitutive  equations  is  that  associated  with 
plane  isotropy. 

Evaluations  of  the  material  are  usually  performed  in  two  directions  normally 
identified  as  "with  the  grain"  and  "against"  or  "across  the  grain.  "  It  has  been 
stated  that  seven  interrelated  elastic  constants  (five  independent  constants) 
describe  the  stress -strain  behavior  of  the  graphite (9, 10),  namely,  three 
different  Poisson's  ratio  {^ee  Figure  10),  two  uniaxial  moduli,  and  two 
shear  moduli. 

At  this  point,  it  is  well  to  emphasize  that  this  mathematical  model  for 
stress -strain  response  has  tne  inherent  limitations  of: 

a.  linear  stress -strain  behavior 

b.  no  distinction  between  tensile  and  compressive  properties. 

In  reality  (see  Reference  9,  page  7),  for  both  grain  orientations,  the  tensile 
and  compressive  properties  differ;  the  compressive  moduli  being  lower  than 
the  tensile  moduli  in  room  temperature  tests  and  higher  in  elevated  temperature 
tests.  As  far  as  lateral  strain  for  uniaxial  loading  is  concerned,  two  of  the 
t  ree  Poisson's  ratios  are  about  the  same  over  a  wide  temperature  range  and 
the  third  is  considerably  higher. 

The  transformation  equations,  which  predict  elastic  constant  variations 
as  a  function  of  orientation  to  the  directions  of  anisotropy,  are  readily 
available  and  stated  in  the  form: 

1  sin4  0  /  1  2  *T2 

E(0)  E2  \G12  E1 
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).  2  2  Cos4  0 

sin  0  cos  0  +  - 

E1 


SECTION  IV. 


EXPERIMENTAL  WORK  AND  DISCUSSION 

A.  Program  Plan 

During  the  past  year  the  program  plan  was  one  of  essentially  two  phases.  The 
first  phase  involved  the  extension  of  the  NDT /mechanical  properties  correlations, 
which  began  during  the  first  year,  to  all  graphite  grades  of  interest  to  the  Air 
Force.  Grades  worked  on  were  Union  Carbide's  ATJ,  RVA,  CFW,  CFZ,  ZTA, 
and  a  special  grade  designated  as  RVA  (-1);  the  RVA  (-1)  billet  was  so  designated 
because  it  was  processed  to  within  one  step  of  completely  processed  RVA.  The 
second  phase  of  the  program  was  concerned  with  developing  correlations  between 
NDT  and  other  properties  and  behavior  characteristics  for  these  same  grades. 

As  in  the  past,  the  program  plan  called  for  a  three-stage  evaluation  approach; 
the  graphite  billet,  selected  sections  and  specimen  blanks,  and  the  final 
destructive  test  specimen.  Of  course,  actual  correlation  of  destructively 
determined  properties  with  NDT  response  provided  the  experimental  verification 
of  the  applicability  of  the  various  NDT  techniques.  Practical  NDT  measurement 
and/or  inspection  systems  for  implementing  these  correlations  on  graphite  billets 
and  hardware  were  also  a  primary  consideration. 

B.  NDT  Screening 

Screening  of  billets  for  sectioning  and  selection  of  specimen  blanks  to  embrace 
the  maximum  range  of  material  variability  was  performed  by  ultrasonic  velocity 
and  radiometry  measurements.  Both  of  these  NDT  measurements  are  sensitive 
to  variations  in  density,  which  has  been  identified  and  experimentally  verified 
to  be  the  predominant  material  variable  affecting  the  properties  of  interest. 
Radiometry  can  be  used  to  obtain  a  direct  measurement  of  density  to  within  ±1%. 
Velocity  measurements,  in  addition  to  their  sensitivity  to  density,  are  sensitive 
to  the  material  anisotropy  and  to  elastic  properties.  A  combination  of  radiometrically 
determined  density  and  longitudinal  wave  velocity  is  uniquely  correlatable  to  other 
properties  of  interest,  as  will  be  discussed  later.  These  reasons  provide  sufficient 
justification  for  the  use  of  these  readily  applied  state-of-the-art  techniques  for 
the  required  screening.  Moreover,  at  the  present  state  of  development,  other 
NDT  techniques,  such  as  infrared. and  eddy  current,  do  not  supply  significant 
screening  capabilities  beyond  those  stated  above.  For  a  more  detailed  description 
of  screening  techniques,  as  well  as  examples  of  typical  variations  within  and  between 
billets,  the  reader  is  referred  to  last  year's  final  report  (Reference  1). 

C.  NDT /Mechanical  Properties  Correlations 

The  relationship  of  the  longitudinal -wave  velocity  to  dynamic  modulus  and  density 
has  been  discussed  in  detail  in  previous  reports,  and  may  be  expressed  as; 

^(1 -a) 

VL  =  -  (51 

L  (f(l  +<7)(1  -2a) 

where  (Ed).  (  p  land  (  a  )  are  the  dynamic  modulus,  density  and  Poisson's  ratio, 
respectively. 


An  attempt  was  made  to  determine  the  effect  of  Poisson's  ratio  in  this  expression. 
To  do  this  it  was  first  necessary  to  obtain  the  range  of  values  to  be  expected  for  <r 
in  representative  grades.  Poisson's  ratio  was  calculated  on  18  selected  specimens 
of  ATJ  and  10  selected  specimens  of  RVA  from  measurements  of  transverse-wave 
velocity  (  Vy  )  and  longitudinal-wave  velocity  (  VL  )  using  the  relationship: 


1- 2(Vt/Vl)2 

2- 2(Vt/Vl)2 


(6) 


It  should  be  noted  that  the  specimens  selected  were  designed  to  cover  the  widest 
range  of  material  variability  available  in  these  grades,  both  with  and  against  the 
grain.  Although  the  above  expression  is  only  strictly  valid  for  an  isotropic 
material,  the  calculations  obtained  were  felt  to  be  sufficiently  accurate  to 
identify  the  range  of  (  o  )  to  be  expected,  which  in  all  cases  was  calculated  to 
be  less  than  0.  20.  As  added  support,  6  determinations  of  (  a  )  were  made  on 
CFW  graphite  by  compression  testing  (4  against  the  grain  and  2  with  the  grain). 
The  maximum  value  obtained  against  the  grain  was  0. 15,  while  the  maximum 
obtained  with  the  grain  was  0.  11,  which  is  within  the  range  identified  from  the 
calculations  above. 


The  function  1(a)  =  - - -  in  equation  (5)  has  been  graphed  versus  (a) 

(l+o)(l-2o) 

in  Figure  1.  Note  that  for  positive  values  of  (a),  the  function  f  (a)  has  the  limiting 
value  of  1. 00  for  a  =  zero,  and  becomes  infinite  at  a  =  0.  50.  Also,  for  values  of 
(a)  less  than  0.20,  f  (  a)  is  less  than  1.10.  For  our  purposes  then,  equation  5 
may  be  written  as: 


This  expression  may  in  turn  be  written  as: 


log  VL  =  1/2  log  Ed  -  1/2  log  p  ^ 

For  purposes  of  analyzing  the  data,  it  is  sometimes  more  convenient  to  work 
with  equation  8  since  log -log  displays  of  the  data  involve  straight  line  relationships. 
Figure  2  is  a  graph  of  such  a  plot  for  Union  Carbide  grades  ATJ,  RVA,  CFW, 
and  CFZ.  The  data  from  any  given  billet,  with  and/or  against  the  grain,  fits  a 
given  line;  the  ATJ  data  for  both  directions,  by  virtue  of  its  lower  density,  on 
the  upper  line  and  the  other  grades  of  higher  density  on  the  lower  line.  The 
sensitivity  to  approximately  a  5  percent  difference  in  average  density  between 
ATJ  and  the  other  grades  is  noted  to  be  a  significant  vertical  displacement  in 
the  graphs.  Furthermore,  the  slopes  of  the  linei?  are  approximately  1/2  as 
predicted  by  equation  (8). 
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Fijur.  2  LONGITUDINAL  VELOCITY  (VL)  VERSUS  DYNAMIC  MODULUS  <Ej) 
FOR  GRADES  ATJ,  RVA,  CFW,  AND  CFZ 


Another  useful  way  of  presenting  this  data  is  that  shown  in  Figure  3.  Here, 

U' V| )  is  plotted  as  a  function  of  the  dynamic  modulus  (Ed)  and  is  a  more  useful 
graph  since  nondestructively  measurable  quantities  (  p  and  V^)  are  plotted 
versus  a  property  of  interest.  Figures  2  and  3  provide  the  experimental 
verification  of  the  relationship  shown  in  equations  (7)  and  (8).  The  log-log 
display  was  presented  to  illustrate  the  effect  of  density  since  it  has  been 
identified  as  the  first-order  material  variable  which  influences  properties 
and  behavior  characteristics  of  interest.  Practical  aspects  are  that  the 
independent  nondestructive  measurement  of  velocity  and  density  can  be  used 
to  determine  dynamic  elastic  modulus,  both  with  and  against  the  grain  to 
within  approximately  ±0.06  x  10^  psi  in  billets  and  hardware  components 
of  these  grades  of  graphite. 

The  data  from  grades  ZTA  and  RVA  (-1)  were  deliberately  omitted  from 
the  graphs  of  Figures  3  and  4  because  the  range  of  values  far  exceeded  those 
of  the  other  grades  plotted.  For  the  16  ZTA  specimens  shipped  to  Avco 
for  evaluation  (5  specimens  with  the  grain  and  10  against)  values  of  dynamic 
modulus  ranged  from  0.  77  x  10^  psi  against  the  grain  to  3.  39  x  10^  psi  for 
specimens  with  the  grain.  Specimens  of  RVA  (-1)  (2  with  the  grain  and  2  against 
the  grain):  ranged  from  1.  83  x  10^  psi  to  2.  71  x  10^  psi.  To  illustrate  the  range 
of  these  values  compared  to  the  other  grades  tested  and  plotted  in  Figure  3, 
consider  the  graph  of  Figure  4.  The  solid  line  is  the  graph  of  the  data  for  the 
grades  shown  in  Figure  3,  while  the  dotted  line  represents  the  extension  of  the 
graph,  with  and  against  the  grain,  to  include  the  ZTA  data.  It  is  significant 
to  note  that,  within  the  same  scatter  as  shown  for  the  data  in  Figure  3,  the 
ZTA  data  fits  the  same  line.  The  change  in  intercept  for  the  RVA  (-1)  data 
has  not  been  explained.  However,  it  has  been  noted  that  dynamic  modulus 
measurements  by  the  resonance  beam  technique  were  difficult  to  obtain  on 
this  material  because  of  poor  definition  of  the  resonant  frequency. 

Having  completed  all  the  required  nondestructive  testing  on  the  specimen 
blanks  for  the  various  billets  of  each  grade,  selected  specimen  blanks  covering 
the  desired  range  of  NDT  variability  were  machined  into  tensile  bars  for 
destructive  mechanical  evaluation.  Longitudinal-wave  velocity  was  measured 
along  the  axis  of  each  tensile  specimen  to  obtain  data  in  the  direction  of  loading 
prior  to  tensile  testing.  Radiographs  of  each  tensile  specimen  were  also  taken 
to  assess  macroscopic  flaw  content  and  thereby  avoid  extraneous  data  which 
might  be  caused  by  flaws. 

Establishing  an  NDT/tensile  strength  correlation  involved  some  rather  unique 
considerations.  To  illustrate  these  ideas,  consider  Figure  5,  which  shows  a 
stress-strain  curve  for  a  typical  specimen  from  a  given  direction  in  a  typical 
billet  of  graphite.  Also  shown  in  the  figure  is  the  graph  of  a  line  whose  slope  is 
numerically  equal  to  the  product  (pVL)  for  that  specimen.  The  value  of  density 
in  tnis  case  was  taken  from  a  section  of  the  tensile  specimen  adjacent  to  the 
fracture  and,  as  noted  earlier.  (VL)  was  measured  along  the  axis  c^f  the  specimen 
prior  to  loading.  It  will  be  remembered  that  for  our  purposes  (pVL)  was  noted 
to  be  equal  to  the  dynamic  modulus  (Eq).  Figure  5  also  shows  the  secant  modulus 
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Figur*  5  TYPICAL  STRESS-STRAIN  PLOT  SHOWING  RELATIVE  POSITIONS  OF  (pV2  ), 
TENSILE  MODULUS  (Ej)  AND  SECANT  MODULUS  (E,)  FOR  A  TYPICAL 
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taken  from  origin  to  failure  on  the  stress -strain  curves.  The  secant  modulus 
is  defined  here  as*. 

Ultimate  Tensile  Stress  (UTS) 

*  Total  Stain  to  Failure  (r*f ) 

Analysis  of  the  tensile  data  revealed  that  the  ratio  of  (pV^)  to  (  F-s  )  was 
essentially  a  constant  for  a  given  direction  in  a  given  grade  of  graphite.  Figure  6 
is  a  plot  of  (p Vj2 }  versus  the  secant  modulus  {  Es  )  for  specimens  taken  with  and 
against  the  grain  from  grades  ATJ,  RVA,  CFW  and  CFZ.  Actual  ultimate  tensile 
stress  and  total  strain-to-failure  on  each  individual  specimen  were  used  to 
calculate  (Es).  Since  graphite  is  brittle  at  room  temperature,  variations  in 
total  strain  between  specimens  of  the  same  direction  in  a  given  grade  are  noted 
to  be  small.  Therefore,  for  purposes  of  determining  ultimate  tensile  strength 
from  the  graph,  the  total  strain  for  specimens  from  the  same  direction  of  the 
same  grade  is  taken  to  be  constant;  the  average  values  being  used  in  each  case. 

In  practice,  the  density  (  p  )  and  the  velocity  (VL)  are  nondestructively 
measurable  quantities  which  can  be  obtained  on  billets  and  hardware.  A  given 
value  of  (pv£)  corresponds,  within  the  limits  shown,  to  the  secant  modulus.  By 
multiplying  the  secant  modulus  by  the  average  total  strain  for  a  given  direction 
of  a  given  grade  of  graphite,  the  ultimate  tensile  strength  may  be  determined  to 
within  approximately  ±300  psi.  This  range  is  well  within  the  range  of  values 
normally  reported  for  this  property  in  these  grades  at  room  temperature. 

Again,  the  data  from  grades  ZTA  and  RVA  (-1)  was  deliberately  omitted 
from  the  graph  of  Figure  6  because,  as  noted  for  the  pv(  ^/dynamic  modulus 
correlation  of  Figure  4,  the  range  of  valuer  far  exceeded  that  of  the  other 
grades.  Secant  modulus  ranged  from  as  low  as  ab'-ut  0.  70  x  10^  psi  against 
the  grain  to  about  3.  4  x  10^  psi  with  the  grain,  which  is  approximately  3  times 
the  range  noted  for  the  other  grades  tested.  Specimens  of  RVA  (-1)  ranged 
from  1.  25  x  10^  psi  against  the  grain  to  1.98  x  10^  psi  with  the  grain.  Figure 
7  shows  the  values  of  data  for  the  grades  plotted  i Figure  6,  while  the  dotted 
line  represents  the  extension  of  the  correlation,  with  and  against  the  grain, 
to  include  the  ZTA  and  RVA  (-1)  data.  It  is  significant  to  note  that,  as  was 
the  similar  case  for  the  dynamic  modulus  correlation,  the  ZTA  and  RVA  (-1) 
data  fits  the  same  line,  which  supports  the  validity  of  the  analysis  which  led 
to  these  correlations. 
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D.  Preliminary  Investigation  of  the  Thermo-Mechanical  Behavior  of  Graphite 

(7  8) 

Several  investigators  ’  have  suggested  that,  because  of  the  complex 
combinations  of  failure  mechanisms  that  may  occur,  prediction  of  the 
behavior  of  potential  rocket  nozzle  materials  based  on  physical  property 
data  or  simple  laboratory  tests  is  difficult,  if  not  impossible. 

For  graphite,  this  is  in  part  due  to  ignorance  of  the  true  failure  mechanisms 
of  such  materials  and  also  partly  attributable  to  a  general  over-confidence 
in  simple  mathematical  idealization  5f  graphitic  mechanical  and  thermal 
behavior.  In  this  discussion,  the  widely  applied  assumption  of  plane  isotropy 
is  first  presented,  and  an  effort  is  then  made  to  delineate  some  of  the 
controlling  influences  of  actual  response  phenomena. 

Typicai  microstructure  of  molded  ATJ  graphite  is  illustrated  by  the 
photomicrographs,  Figures  8  and  9.  The  non-uniform  density  as  well  as 
the  inhere:*;  porosity  is  clearly  revealed  by  these  magnified  views  of  the 
material.  Due  to  observed  anisotropy  and  the  nature  of  the  molding  process, 
a  widely  suggested  form  of  constitutive  equations  is  that  associated  with 
plane  isotropy. 

Evaluations  of  the  material  are  usually  performed  in  two  directions  normally 
identified  as  "with  the  grain"  and  "against"  or  "across  the  grain.  "  It  has  been 
stated  that  seven  interrelated  elastic  constants  (five  independent  constants) 
describe  the  stress -strain  behavior  of  the  graphite'^*  10) t  namely,  three 
different  Poisson's  ratio  {i?ee  Figure  10),  two  uniaxial  moduli,  and  two 
shear  moduli. 

At  this  point,  it  is  well  to  emphasize  that  this  mathematical  model  for 
stress -strain  response  has  tne  inherent  limitations  of: 

a.  linear  stress-strain  behavior 

b.  no  distinction  between  tensile  and  compressive  properties. 

In  reality  (see  Reference  9,  page  7),  for  both  grain  orientations,  the  tensile 
and  compressive  properties  differ;  the  compressive  moduli  being  lower  than 
the  tensile  moduli  in  room  temperature  tests  and  higher  in  elevated  temperature 
tests.  As  far  as  lateral  strain  for  uniaxial  loading  is  concerned,  two  of  the 
t  ree  Poisson's  ratios  are  about  the  same  over  a  wide  temperature  range  and 
the  third  is  considerably  higher. 

The  transformation  equations,  which  predict  elastic  constant  variations 
as  a  function  of  orientation  to  the  directions  of  anisotropy,  are  readily 
available  and  stated  in  the  form: 


sin^  cf> 


2  v 


12 


:<0) 


4  . 

.  2  2  Cos  <p 

sin  (f>  cos  <f>  +  — - 


M2 


-1 


(10) 
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Figur.9  TYPICAL  MOLDED  ATJ  GRAPHITE  STRUCTURE,  120X,  X  Z  PLANE 
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Figure  !0  ILLUSTRATION  OF  THREE  POISSON'S  RATIO  IN  A  MOLDED  GRAPHITE 


where:E(<£)is  the  uniaxial  modulus  at  angle  <j>  to  the  Xj  direction, 

Ej  ,  E2  are  the  uniaxial  moduli  in  X! ,  X2  directions 
G12  is  the  appropriate  shear  modulus 
t>12is  the  Poisson's  ratio 

1  /  1  1  2  v\2  1  \  ,  ,  1 
— — —  =  4  1—  +  —  +  —  - - -  I  sin*  4>  cos*  <f>  +  .  -  ■ 

G12(<£)  \E1  E2  E1  G12  /  G12 

However,  a  relatively  thorough  literature  survey  did  not  yield  any  information 

as  to  the  validity  of  such  a  predicted  properties  variation.  Preliminary 

studies  were  completed  in  this  area.  Eight  torsion  samples  were  machined 

from  a  typical  ATJ  billet,  and  velocity  measurements  made  as  listed  in  Table  1, 

with  the  associated  dynamic  moduli  determined  from  correlation  data  reported 

in  Section  IV-C. 

TABLE  I 


Uniaxial  Modulus  Properties  of  the  Torsion  Samples 


Sample 

Velocity 

Number 

Direction 

in Ifi  sec. 

Approximate  Dynamic  Modulus,  psi 

1 

Y 

.1047 

1.68  x  106 

2 

Y 

.1058 

1.19  x  106 

3 

X 

.0888 

4 

X 

.0886 

/ 

5 

Z 

.1009 

1.  52  x  10° 

6 

X 

.1017 

/ 

7 

45 e  in  YZ 

.1020 

1.58  x  10° 

8 

plane 

.1020 

Inspection  of  the  dynamic  modulus  values  suggests  that  the  YZ  plane  departs 
from  assumed  isotropy. 

The  torsion  test  data  is  listed  in  Table  II.  Sample  geometry  and  test 
fixtures  are  shown  in  Figure  11  and  typical  torsion  failures  illustrated  by 
Figure  12. 

TABLE  II 

Torsion  Test  Results 


Sample  No. 

Axis 

Shear  Modulus 

Shear  Strength 

1 

Y 

6 

.  641  x  10  psi 

3830  psi 

2 

Y 

.678  x  106 

3960 

3 

X 

. 569  x  106 

3560 

4 

X 

.  561  x  106 

3750 
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TABLE  II  (Cont'd) 


Sample  No.  Axis 


Shear  Modulus 


Shear  Strength 


5 

6 

7 

8 


Z 

Z 

45®  in  YZ 
45°  in  YZ 


.  625  x  10^  psi 
.  613  x  106 
.  643  x  10* 

.  595  x  106 


3460  psi 
3880 
4070 
3400 


It  has  been  suggested^)  that  torsion  tests  are  not  appropriate  for  shear 
properties  determinations  since  tensile  stresses  f.nduce  the  actual  fracture. 
From  this  point  of  view,  it  is  interesting  to  note  that  the  maximum  shear 
strengths  reported  in  Table  II  are  in  rather  close  agreement  with  data  presented 
in  Reference  9.  In  the  same  reference,  a  notched  bar  shear  test  is  advocated 
for  both  shear  moduli  and  shear  strength  observations.  However,  data  reported 
in  that  study  indicates  gage  length  effects  on  both  modulus  and  shear  strength. 


Therefore,  it  is  currently  b  elieved  that  the  shear  moduli  determined  from 
torsion  tests  are  reliable.  In  order  to  substantiate  the  shear  modulus 
measurement,  two  torsional  resonance  vibration  tests  were  performed  on 
circular  bars  and  the  resulting  shear  modulus  values  were  in  excellent 
agreement  with  the  torsion  test  data. 


Figure  13  illustrates  theoretical  uniaxial  modulus  variation  and  Figure  14 
presents  a  comparison  of  predicted  and  observed  shear  modulus.  Fair 
agreement  seems  evident.  Keeping  in  mind  the  Known  thermal  effects  on 
mechanical  properties  of  the  graphite,  it  is  suggested  that  a  valuable 
experimental  program  would  involve  a  comparison  of  mechanical  response 
at  various  orientations  over  a  wide  temperature  range. 


It  has  been  demonstrated  by  means  of  wave  velocity  measurements  that 
departures  from  isotropy  exist  in  the  YZ  plane,  that  the  shear  modulus  variation 
approximates  theoretical  predictions. 


One  of  the  goals  of  this  research,  as  mentioned  earlier,  was  to  obtain  an 
empirical  NDT  parameter  for  assessing  the  resistance  of  graphite  material  to 
thermal  shock.  It  has  been  demonstrated  that  the  graphite  behaves  as  an 
anisotropic  medium.  Therefore,  in  order  to  design  meaningful  experiments 
to  formulate  a  more  rational  theoretical  basis  for  the  associated  NDT  measurements, 
thermoelastic  stress  wave  phenomena  must  be  understood.  A  detailed  presentation 
of  the  equation  of  wave  propagation  in  transversely  isotropic  media  appears  in 
the  appendix.  This  full  set  of  field  equations  must  be  resorted  to  in  order  to 
determine  the  importance  of  coupled  thermo-mechanical  behavior. 


For  most  isotropic  elastic  material,  the  coupling  effects  in  the  thermoelastic 
field  equations  are  justifiably  negligible.  However,  for  graphites,  the  order  of 
magnitude  of  the  coefficients  of  the  equations  has  not  been  completely  defined 
over  a  wide  temperatu.e  range  and  so  the  importance  of  the  coupling  terms 
is  somewhat  unsettled. 
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Beyond  these  matters,  complexities  associated  with  heat  conduction  in 
which  the  body  undergoes  a  phase  change  during  exposure  to  the  thermal 
environment  undoubtedly  are  important  in  relation  to  thermal  shock  phenomena. 
These  reasons  all  add  credence  to  the  statement  that  full-scale  performance 
of  rocket  nozzles  cannot  be  predicted  from  ordinary  properties  data. 

However,  this  is  only  in  part  true  since  the  refinements  of  more 
sophisticated  mathematical  representations  have  not  been  applied. 

Admittedly,  even  the  asymptotic  solutions  to  these  problems  are  difficult, 
time  consuming  and  costly.  Rocket  nozzle  failures,  however,  are  also 
costly  and  dangerous  as  well. 

T  vs  failure  response  of  isotropic  materials  can  be  considered  from  both 
the  microscopic  and  macroscopic  point  of  view.  From  the  viewpoint  of 
gross  engineering  types  of  response,  it  is  useful  to  enumerate  the  ordinarily 
used  fracture  criteria  for  elastic  media. 

1.  maximum  principal  stress 

2.  maximum  principal  strain 

3.  maximum  principal  stress  difference  (shear  stress) 

4.  maximum  principal  strain  difference  (shear  strain) 

5.  maximum  distortional  strain  energy 

6.  maximum  conserved  distortional  strain  energy 

7.  maximum  total  strain  energy 

The  important  point  to  note  is  that  no  universal  fracture  criterion  has  been 
established  and  that  the  success  of  a  given  fracture  hypothesis  depends  in  large 
measure  upon  the  material  with  which  it  is  associated. 

One  can  quickly  list  the  various  yield  criteria  and  attempt  to  verify  the 
applicability  of  these  theories  to  graphite  behavior.  For  example,  the  VonMises 
yield  criterion  based  on  the  mean  value  of  the  principal  stress  difference  takes 
the  form 


"  3 )  t-  (aj-wjT 


(U) 


for  simple  tension  this  equation  results  in  oQ  ,  while  for  pure  shear,  the  yield 
stress  =  (’0/\T.  In  graphites,  the  shear  strength  is  of  the  same  magnitude  as 
the  tensile  strength  which  negates  such  simple  criterion. 

Hill^^  postulated  that  the  yield  condition  for  anisotropic  media  is  a 
quadratic  function  of  the  stress  components 

2f(ojj)  -  +■  f  H(f/j  ct2)2 

f  2  L  r2f  t  2Mr32j  +  2Nr12  ^ 

where  Xj,  X2 ,  X3  are  axes  of  material  symmetry  and  F,  G  ,  H  ,  L  ,  M  ,  and  N 
are  material  constants. 
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The  strength  of  quasi -homogeneous  anisotropic  composites  was  reported 
by  Azzi  and  Tsai^^'.  One  of  the  basic  assumptions  of  their  failure  condition 
is  that  there  exist  three  mutually  perpendicular  planes  of  symmetry  within  the 
anisotropic  body  and  the  material  is,  therefore,  taken  to  be  orthotropic 
strengthwise.  Under  this  assumption  the  coordinate  system  for  the  state  of 
stress  and  the  material  symmetry  must  be  the  same.  Thus,  the  state  of 
stress  imposed  on  a  body  must  be  transformed  to  the  coordinate  system  of 
material  symmetry  and  then  the  yield  conditions  applied.  For  uniaxial  ultimate 
strength,  Tsai  has  suggested  that  the  applicable  equation  is 


The  appropriateness  of  such  an  equation  for  graphite  materials  has  not 
been  investigated.  However,  an  alternate  equation  may  be  derived  for  the 
case  of  pure  shear  (see  Reference  13,  page  12).  For  the  case  where 
<f>  =  45°  ,  the  shear  strength  takes  the  form 


(14) 


which  reduces  to  o0/\J~T  when  ctj  =  ct2  and  therefore,  agrees  with  the  Von  Mises 
condition.  The  applicability  of  this  equation  for  predicting  shear  strength  can 
be  assessed  by  using  the  data  of  Table  II  as  well  as  that  available  in  the 

literature^. 


Considering  the  lack  of  success  of  simple  yield  criterion,  and  the  fact 
that  correlations  with  nondestructive  tests  can  be  established  between 
modulus,  strength  and  porosity  of  graphite,  it  is  suggested  that  a  more 
realistic  yield  criterion  might  be  the  Mohr  Coulomb  theory  ordinarily 
applied  to  granular  materials.  This  is  particularly  true  since  in  small 
solid  propellant  rocket  engines  applications  the  nozzle  insert  is  in  part 
subjected  to  a  hydrostatic  compression. 

The  failure  condition  can  be  stated  in  the  form 


rmax 


k  +  an  tan  <f> 


(15) 


and  the  proof  of  such  behavior  can  be  established  by  performing  triaxial 
tests  on  cylinders  of  the  material.  For  this  purpose,  the  apparatus 
illustrated  in  Figure  15  was  used  to  subject  membrane  covered  molded 
ATJ  graphite  to  several  levels  of  hydrostatic  pressure  and  then  to  compress 
the  sample  to  failure  with  the  axial  loading  piston.  A  detailed  discussion 
of  routine  and  special  triaxial  test  techniques  appears  in  Reference  14.  For  the 
ATJ  specimens,  the  experiments  were  performed  on  both  with  and  against  the 
grain  orientations  and  the  results  are  presented  in  Figures  16  and  17. 
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These  yield  envelopes  show  that  the  with  the  grain  orientation  has  significant 
hydrostatic  pressure  effects  on  the  maximum  shear  strength  while  only  very 
slight  influences  were  evident  for  the  against  the  grain  orientation.  This 
suggests  that  the  with  grain  orientation  exhibits  response  characteristics 
similar  to  granular  media. 

In  order  to  evaluate  the  influence  of  this  yield  phenomena  a  simple  plasticity 
problem  is  now  considered  by  comparison  of  elementary  and  Mohr  Coulomb 
criteria  for  the  case  of  a  rigid  plastic  wedge.  The  solutions  for  this  case 
with  uniform  normal  pressure  on  one  face  are  compared  for  each,  (see  Figure  18) 

Reference  15  lists  the  solution  for  the  simple  yield  criteria  as 


P  =  2k  (1  +  2  /80  -  90°) 


(16) 


,(16) 


while  for  the  case  of  the  Mohr  Coulomb  Criteria,  Sokolovsky'  '  suggests 
substituting  the  equations 


W  =  k  +  °n  tan  <t> 


°\-a2 


°n  =  y  (al  +  *2>  -  l - - 


sin  c6 


(17) 


into  the  yield  criterion  F  (<7”)  so  that  one  has 

T  1  1  2  1  „  .  2  1 

k  cos  +  —  (oj  +  o^)  sin  ^  —  (<7j  -  Oj)  cos  cf>  -  —  (rq  -  nj)  sin  <pl 


(18) 


Satisfying  the  equilibrium  equations  then  leads  to  the  expression  for  the 
limiting  plastic  pressure  on  the  wedge 


p '  =  k  cot  <f> 


2  ft,  >  7 


■,  /  IT  <f>\ 

exp  (j  2/3q  -  rr  I  tan  cf>)  tan*  ^4”  +  ~2J~ 


(19) 


The  ratio  of  p'/p  determined  using  equations  16  and  17  where  k  and  <j>  are 
obtained  from  Figure  16,  are  compared  in  Figure  19  as  a  function  of  wedge 
angle.  The  solution  of  course  is  not  valid  for  angles  approaching  180°  since 
this  geometry  would  lead  to  a  problem  of  contained  plasticity.  In  the  range  of  90 ‘ 
to  140°  ,  the  observed  ratio  shows  a  variation  of  20  to  30%, 

It  is  suggested  that  the  Mohr  Coulomb  Criteria  is  mathematically  tractable 
and  is  perhaps  a  more  realistic  description  of  failure  of  the  graphite  under 
ce  rtain  conditions . 


-33- 


E.  Investigation  of  Infrared  Techniques  for  Nondestructive  Determination 
of  Thermal  Conductivity 


The  thermal  conductivity  of  materials  is  an  important  physical  quantity  to 
the  engineer  involved  in  the  design  of  fabrications  requiring  heat  transfer 
consideration.  Conductivity  is  a  basic  quantity  that  is  absolutely  necessary 
for  determining  the  amount  of  heat  that  can  flow  in  a  material  continuously. 

When  combined  with  values  of  the  density  and  heat  capacity  of  a  material 
the  conductivity  is  proportional  to  a  thermal  parameter  called  diffusivity, 
t  y  which  the  spatial  intensity  and  temporal  distribution  of  thermal  transients 
in  the  material  can  be  determined.  When  conductivity  is  combined  with 
knowledge  of  a  material's  elastic  modulus,  tensile  strength  and  expansion 
coefficient,  the  relative  ability  of  high -temperature  refractories,  such  as 
graphite,  to  withstand  thermal  shock  can  also  be  determined. 

Determination  of  a  material's  thermal  conductivity  has  been  traditionally 
performed  under  steady-state  conditions  of  heat  flow  on  specimens  having 
some  experimentally  suitable  geometry.  Some  recent  studies  have  described 
the  use  of  a  transient  heating  through  transmission  technique  to  determine 
conductivity,  diffusivity  and  heat  capacity  _  These  studies  also 

have  required  the  use  of  special  specimen  geometries  and  the  need  for 
access  to  two  opposite  surfaces,  both  of  which  necessitate  the  destructive 
testing  of  a  fabrication.  These  technique  requirements  result  principally 
from  limitations  imposed  in  attempting  to  experimentally  achieve  boundary 
conditions  suitable  for  analysis. 

Indirect  nondestructive  testing  techniques  for  determining  thermal 
conductivity  in  many  electrical  conductors  are  feasible  by  correlating  their 
electrical  and  thermal  conductivities.  Electrical  measurements,  however, 
can  usually  be  performed  only  on  large  sections  and  result  in  bulk  determinations. 
For  other  materials,  such  as  dielectrics  and  graphite,  suitable  correlations 
between  their  thermal  conductivities  and  other  physical  quantities  have  not 
been  established  as  yet.  Nondestructive  direct  measurements,  then,  must 
be  performed  especially  on  these  latter  materials,  and  preferably  in  a  practical 
manner  that  yields  results  in  a  comparatively  localized  region,  in  order  to 
observe  maximum  material  variability  and  contributing  processes. 

It  is  believed  that  prior  to  this  study  the  direct  nondestructive  determination 
of  the  thermal  conductivity  of  fabrications  has  not  been  performed.  To  this 
end,  a  novel  technique  has  been  developed  from  which  the  near-room-temperature 
conductivity  of  a  solid  can  be  determined  directly.  Theory  is  presented  and  is 
supported  by  limited  experimentation  to  indicate  the  feasibility  of  the  technique. 

The  technique  is  simple  and  consists  of  simultaneously  heating  a  surface  of  a 

"semi-infinite"  solid  by  radiation  and  observing  the  temperature  history  at 

this  surface  using  an  infrared  radiometer  as  the  radiance  sensor.  Analysis 

of  the  history  record  by  comparison  with  a  standard  record  and  knowledge  of 

the  unknown's  density  and  specific  heat  permits  the  determination  of  the  conductivity, 
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Type  ATJ  graphite  was  the  material  selected  for  emphasis  in  this  study 
because  of  Avco's  present  interest  in  it,  steming  from  its  range  of  properties 
variability.  There  is  every  reason  to  believe,  however,  that  the  technique 
can  be  extended  to  include  testing  of  solids  in  general. 

1 .  Theory 

The  following  theoretical  development  is  aimed  at  obtaining  the  thermal 
conductivity  of  a'solid  by  obser-dng  the  temperature  history  of  an  exposed 
surface,  using  an  infrared  radiometer,  while  this  surface  is  subjected  to  a 
covstant  radiant  heat  flux.  By  comparing  geometrical  characteristics  of 
this  history  with  that  for  a  solid  having  a  known  conductivity,  together  with 
knowledge  of  the  unknown's  density  and  heat  capacity,  the  unknown  conductivity 
may  be  determined.  This  transient  technique  does  not  require  knowledge 
of  the  flux  incident  to  the  surface,  which  is  a  most  desireable  consequence. 

For  validity,  however,  it  does  require  adherence  to  certain  boundary  conditions 
and  to  criteria  for  nondestructive  testing  and  observability  of  the  history  recording. 

2  .  Experimental  Parameters 


a  .  Fundamentals  of  Heating  Method 

Three  criteria  are  used  as  bases  in  this  development.  The  first 
criterion  presumes  the  solid  to  be  dimensionally  large  enough  so  that 
it  appears  semi-infirite  to  heat  flow  during  the  period  of  simultaneous 
heating-observation.  If  it  is  assumed  that  the  heat  flow  front  is  spatially 
uniform  in  intensity  ar.d  is  unidirectional,  so  that  no  lateral  heat  flow 
loss  occurs,  this  criterion  is  satisfied  for  the  thickness  dimension  by 
the  condition^ 


>  (3at)1/2  . 


(20) 


where  L,*,  is  the  distance  into  the  solid  beyond  which  semi-infiniteness 
exists,  a  is  thermal  diffusivity,  and  t  is  time. 

The  second  criterion  presumes  that  the  temperature  history  of  the  ( 

solid's  surface  is  of  parabolic  form  and  is  generated  by  the  following  relationship 


T(Q.t)  -  T0 


(21) 


where  T  is  the  temperature  of  the  surface  above  an  initial  ambient 
temperature  T0,  Q  is  the  heat  flux  per  unit  area  per  unit  time  absorbed 
at  the  surface,  and  k  is  thermal  conductivity.  The  diffusivity  is  defined 
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in  terms  of  k,  the  density  p  ,  and  the  heat  capacity  at  constant  pressure,  C_  ,  as 


The  third  criterion  presumes  that  the  surface  temperature  of  the  solid 
rises  no  more  than  about  10  degrees  Centigrade  during  the  period  of  heating 
so  that  the  modified  Stephart-Boltzmann  law 


w  =  —  <xa(T4-T04), 


which  describes  the  rate  of  heat  loss  from  a  body  at  T  to  a  body  at  T0,  can 
be  approximated  by  the  relationship 


w  *-(T„  Tq  (T  -  Tc)  . 


In  Equations  (23)  and  (24)  W  iB  the  radiant  energy  or  radiance  emitted  from 
the  surface,  assuming  a  cosine  or  Lambert  law  distribution,  q  is  the 
total  radiant  emissivity  of  the  surface,  and  a  is  the  Stephan-Boltzmann 
constant.  These  two  equations  are  plotted  in  Figure  19,  from  which  it  is 
seen  that  a  10  degree  C  rise  in  the  temperature  of  a  surface  initially  at 
27  degree  C  results  in  a  4  percent  difference  between  the  theoretical  and 
approximate  forms  of  the  radiation  laws.  For  most  solids,  this  10  degree  C 
rise  above  ambient  will  caust  no  significant  change  in  their  physical  properties, 
so  that  testing  is  indeed  nondestructive. 

Since  a  radiometer  is  used  to  detect  the  surface  temperature,  the 
radiant  energy  entering  the  radiometer  is  converted  into  an  electrical 
voltage,  V,  such  that 


V  K  W' 


where  K  is  a  constant  that  accounts  for  this  conversion,  and  W^ 
where  T0^  is  the  radiometer's  reference  source  temperature. 
The  difference  between  tJ  and  TQ  is  a  constant  temperature  5 
temperature  history  curve  is  Actually  one  of  voltage  vs.  time. 


—  <T°To  (T-To' 


Thus  the 
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Figure  20  THEORETICAL  VERSUS  APPROXIMATE  FORM  RADIATION  LAWS 
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Lastly,  in  this  basis  of  development,  Q  is  the  flux  of  heat  deposited 
at  the  solid's  surface  from  the  radiation  source  which  emits  a  flux  F 
incident  to  the  surface.  These  two  quantities  are  related  by  the  equation 


Q  -  F  , 


(26) 


where  a-j-  is  the  total  radiant  absorptaace  of  the  surface.  For  opaque 
gray  surfaces,  tj  -  .  Furthermore,  since  heating  of  the  surface 

is  due  principally  to  that  part  of  the  incident  radiation  energy  spectrum 
in  the  infrared  wavelength  region,  as  is  the  energy  spectrum  emitted 
by  the  heated  solid,  then  both  spectra  correspond  to  the  spectral  region 
detected  by  the  radiometer,  so  that  (j  -  >  where  is  the  emissivity 

of  the  radiometer's  detectible  infrared  spectrum. 


Combining  Equations  (21,22  24-26)  and  solving  for  V  yields  the  relationship 


4  \l/2 


V  =  A 


\ypc  c 


tl/2  _  B 


(27a) 


where 


A  = 


8  K  F  cr  Tq'  ' 
„3/2  ” 


(27b) 


B  = 


41WffTo3s 


(27c) 


V'  = 


V  +  B 


(27d) 


Equation (27a)  is  a  voltage  history  of  parabolic  form  corresponding  to  the 
surface  temperature  history.  The  constants  A  and  B  are  independent  of 
any  characteristics  of  the  solid.  The  latus  rectum  of  the  parabola  equals 

A2  eK 

-  ,  so  that  the  shape  of  the  history  curve  is  dependent  on  the  flux,  F, 

kPcP 

the  emissivity,  ,  and  the  quantity  (kpCp).  A  higher  flux,  or  emissivity, 
or  a  lower  conductivity  increases  the  signal  voltage  and  the  width  of  the 
parabola. 
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temperature  through  a  figure  of  merit  of  the  radiometer  called 
NET=  *  The  NET  is  the  temperature,  seen  by  the  radiometer's 
field  of  view,  required  to  obtain  a  I  =1.  NET  has  units  of 
degrees  Centigrade,  where  small  values  indicate  high  radiometer 
sensitivity  and  performance.  The  NET  for  the  radiometer,  as 
used  in  these  studies,  was  0.  5  degree  C,  at  a  reference  temperature 
of  300  degree  K. 


c .  Minimum  Radiation  Flux 


By  combining  Equations{27)  and  (3^  ,  it  is  possible  to  obtain  a  relationship 
for  the  minimum  incident  radiation  flux,  F0,  necessary  to  obtain  an 
observable  radiometer  signal  voltage  within  the  response  time  of  the 
source-detector-recorder  system.  This  relationship  is 


Fft  (cal. /cm2-  sec) 


Values  of  the  parameters  for  the  radiometer  used  in  thia  study  arei 


Vo 

= 

0.  002  Vv  Its  (rms) 

K 

= 

26.7  volts  -  cm2-sec.  /cal. 

To' 

304*  K 

S 

— 

4*  K 

a 

= 

1.  37  x  10-*2  cal.  /cm^-sec. 

A  signal-to-noise  ratio  2  =  3  will  be  assumed.  Graphite  and  copper  are 
two  solids  considered  and  both  are  presumed  to  be  coated,  each  having  an 
emissivity  value,  ,  equal  to  0.85.  The  quantities  (kpC p  )  are  calculated 
from  published  data  for  k  ,  o  ,  and  Cp  ,  which  are  tabulated  in  Table  I,  and 
are  0.82  and  0.081  cal.  ^/cm  -sec-*  C2  for  copper  and  graphite  (in  a  direction 
with  the  grain),  respectively.  Substitution  of  these  values  into  Equation  (35) 
reduces  it  to 


F 


o 


11.9 


(kpC  p)1/2 

t1/2 


(36) 


For  copper  and  graphite  and  for  a  system  response  time  equal  to  8C  msec, 
Equation  (36)  yields 

F0  =  37.7  cal./cmJ-sec 

cu  =  161  watts /cm^ 

F0g  =  11. 8  cal. /cm^-sec 

=  49  watts/ cm^ 


*NET  is  defined  as  thn  temperature  difference  required  to  produce  an  rms 
signal  voltage  equal  to  the  noise  voltage,  Vn  ,  in  a  1  cps  bandwidth  for  ? 
stated  temperature  in  the  field  of  view  of  the  radiometer. 
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It  is  fortunate  that  Cp  is  constant  for  this  case,  since  the  densities  can 
be  readily  determined  using,  for  example,  a  nondestructive  testing  gamma-ray 
radiometric  density -gauging  technique,  or  an  ultrasonic  velocity  technique. 

b.  Observation  Sensitivity 

1.)  Reflected  Radiation 

All  surfaces  reflect  radiation  to  some  degree,  so  that  >  0  .  This 
creates  a  problem  experimentally.  During  heating  a  part  of  the  radiant 
energy  incident  to  the  solid's  surface  is  reflected  into  the  optics  of  the 
radiometer.  The  total  radiation  sensed  by  the  radiometer  is  composed 
of  this  reflected  radiation  and  the  solid's  self-emitted  radiation.  If 
is  independent  of  temperature,  the  reflected  component  is  constant 
during  heating.  This  effect  can  be  stated  mathematically  as 


V"  =  V  +  VR 


(33) 


where  V  is  the  self-emitted  signal  voltage  described  by  Equation^),  and 
VR  is  the  signal  voltage  caused  bv  the  reflected  radiation.  It  is  to  be 
noted  that  tan  9  ,  from  Equation  (28) ,  is  not  affected  by  Vr  since  it 
depends  on  voltage  differences.  To  observe  the  parabolic -shaped 
informational  signal  voltage,  which  rides  on  top  of  Vr  ,  the  voltage 
difference,  V"-  VR,  must  be  greater  than  the  discrimination  level  of 
the  radiometer.  If  F  and/or  ^  is  too  small,  or  k  too  large,  such 
thatV"-VR  cannot  be  observed,  then  the  heating  curve  will  appear  almost 
square  in  shape,  but  with  slightly  canted  skirts.  These  skirts  are  a 
consequence  of  both  the  finite  response  speed  of  the  radiometer's  detector 
and  the  pulse  rise  time  of  the  radiation  incident  to  the  surface.  Of  course, 
if  Vr  can  be  suppressed  electronically  in  the  radiometer  and  sensitivity 
can  be  increased  to  observe  just  V,  reflection  is  not  a  problem. 

2  Signal-to-Noise  Ratio  (Noise  Equivalent  Temperature) 

This  increase  in  sensitivity  may  introduce  another  discrimination 
problem,  which  invo.ves  both  the  noise  equivalent  temperature  (NET)  of 
the  radiometer  and  the  criterion  of  acceptable  si gnal-to -noise  ratio,  2  , 
for  observability.  The  observable  signal  voltage,  VD,  is  related  to  2 
by  the  expression 

v0  =  vN  a2  -d1/2  ,  (34) 

where  2  =  V/VN  and  VN  is  the  noise  voltage.  If  2  =  1,  then  VQ  =  0  and  an 
informational  signal  cannot  be  observed.  Thus  2  must  be  greater  than 
unity  as  one  criterion  for  detectability.  Values  of  2  >  3  indicate  easily 
observable  informational  signals.  Furthermore,  VN  can  be  related  to 


temperature  through  a  figure  of  merit  of  the  radiometer  called 
NE!\  *  The  NET  is  the  temperature,  seen  by  the  radiometer's 
field  of  view,  required  to  obtain  a  2  =  1.  NET  has  units  of 
degrees  Centigrade,  where  small  values  indicate  high  radiometer 
sensitivity  and  performance.  The  NET  for  the  radiometer,  as 
used  in  these  studies,  was  0.  5  degree  C,  at  a  reference  temperature 
of  300  degree  K. 

c  .  Minimum  Radiation  Flux 

By  combining  Equations{27)  and(34  ,  it  is  possible  to  obtain  a  relationship 
for  the  minimum  incident  radiation  flux,  F0,  necessary  to  obtain  an 
observable  radiometer  signal  voltage  within  the  response  time  of  the 
source -detector-recorder  system.  This  relationship  is 

,1/2  (kpCp)1^  nVQ  I 

F  (cal./cm2-sec)  =  — - -  I- - — - +  • 


2  tl/2 


>i/2  r 

7o  r 


KaT0'3 


Values  of  the  parameters  for  the  radiometer  used  in  thia  study  arei 

Vo  =  0.002  v.  Its  (rms) 

K  =  26.  7  volts  -  cm^-sec.  /cal. 

Tc'  =  304°  K 

8  =  4°K 

a  =  1.  37  x  10"*2  cal.  /cm2-sec. -*  K  . 


A  signal-to -noise  ratio  2  =  3  will  be  assumed.  Graphite  and  copper  are 
two  solids  considered  and  both  are  presumed  to  be  coated,  each  having  an 
emissivity  value,  < ^  ,  equal  to  0.85.  The  quantities  (kpCp  )  are  calculated 
from  published  data  for  k  ,  o  ,  and  Cp  ,  which  are  tabulated  in  Table  I,  and 
are  0.  82  and  0.081  cal.  ^/cm  -sec-*  C2  for  copper  and  graphite  (in  a  direction 
with  the  grain),  respectively.  Substitution  of  these  values  into  Equation  (35) 
reduces  it  to 


=  11.9 


(kpCp)1 


For  copper  aad  graphite  and  for  a  system  response  time  equal  to  8C  msec. 
Equation  (36)  yields 

F0  =  37.  7  cal. /cm "-sec 

cu  =  161  watts /cm^ 

F0(,  =  11. 8  cal. /cm^-sec 

=  49  watts/cm^ 


*NET  is  defined  as  tho  temperature  difference  required  to  product  an  rms 
signal  voltage  equal  to  the  noise  voltage,  Vfj,  in  a  1  cps  bandwidth  for  a 
stated  temperature  in  the  field  of  view  of  the  radiometer. 
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That  is  to  say,  L6L  watts/cm  of  infrared  radiation  incident  to  the 
surface  of  copper,  which  has  an  infrared  emissivity  of  0.  85,  is 
necessary  to  raise  the  temperature  of  the  surface  by  an  amount 
sufficient  for  observation  with  a  signal-to-noise  ratio  of  3,  after  a 
period  of  80  msec,  has  elapsed  from  pulse  initiation.  Similarly, 

49  watts /cm*"  are  necessary  when  studying  graphite  oriented  in  a 
direction  with  the  grain.  These  values  for  F0  are  reasonable  when 
compared  with  the  radiation  source  actually  used  in  this  study.  It 
must  be  remarked  that  those  results  have  validity  only  when:  (1)  the 
solid  appears  semi-infinite  in  thickness  during  the  observation 
period;  (2)  the  incident  radiation  is  uniform  in  intensity  over  the 
surface;  (3)  heat  flow  within  the  solid  is  unidirectional  and  perpendicular 
to  the  surface;  and  (4)  the  radiometer's  field-of-view  on  the  surface  is 
located  such  that  the  region  of  observation  experiences  the  three 
foregoing  boundary  conditions.  •  * 

d.  Effect  of  Coatings 

1. )  Emissivity 

The  effect  of  those  coatings  which  produce  an  emissivity  change 
at  the  solid's  surface  on  the  radiation  flux  requirement  can  be  determined 
from  Equation(35)  by  differentiating  it  with  respect  to  and  solving  for 
the  relative  flux  change  AF0/FQ.  This  yields  the  expression: 


The  minus  sign  in  (a)  indicates  that  a  decrease  in  emissivity  requires 
an  increase  in  the  flux  requirement  to  produce  an  observable  signal. 
Substitution  of  the  radiometer  and  S/N  parameters  listed  above  into 
Equation (37b)  gives  the  value  D/E  »  1.5  ,  so  that  Equation  (37b)  reduces  to: 


(<A  +  3)  A,a 
Ux  +  1.5) 


(38) 
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One  example  of  the  use  of  Equation  (38)  is  to  consider  a  5  percent 
decrease  in  «iit  an  =  0.85.  This  emissivity  change  requires  a 
corresponding  increase  in  F0  of  6.7  percent.  Similarly,  a  5  percent 
increase  in  at  =0.1  indicates  that  F0  can  be  diminished  by  9.  7  percent. 

2. )  °ulse  Retardation 


Since  the  solids  considered  in  this  study  are  good  thermal  conductors, 
the  effect  of  a  poor  conducting  coating,  such  as  a  lacquer,  having  a 
thickness  as  small  as  10"^  cm,  is  to  initially  retard  the  flow  of  heat  into 
the  solid(23)4  As  a  consequence  of  this  retardation,  the  surface  temperature 
is  somewhat  greater  during  the  first  few  seconds  after  pulse  initiation 
compared  to  the  temperature  that  would  result  from  an  ideal  coa  ing  having 
the  same  emissivity  but  an  infinite  surface  conductance.  If  the  same  coating 
is  applied  to  solids  having  different  conductivities,  the  rate  of  temperature 
rise  will  be  greater  the  lower  the  conductivity  of  the  solid.  Thus,  the 
tangent  of  the  angle  of  temperature  rise,  in  Equation(28b)  will  be  somewhat 
greater  the  poorer  the  solid  conductor.  For  example,  since  graphite  has 

a  lower  conductivity  than  copper,  the  ratio  of  tangents,  tan  9  (graphite) 

tan  9  (-opper) 

will  be  less  for  a  coating  of  finite  conductance  than  for  the  ideal  case  of 

infinite  surface  conductance.  Similarly,  the  ratio  (kpCp)  (copper) _ 

(kpCp)  (graphite) 

will  be  less  than  for  the  ideal  case.  Moreover,  for  lacquer -coated  graphite, 

the  ratio  _ k  (with  the  grain)  will  also  be  less  than  for  the  ideal  case. 

k  (against  the  grain) 

History  Curve  Analysis 


The  voltage  histories  for  the  several  solids  studied  have  been  calculated 
using  Equation (27)  an  he  published  data  listed  in  Table  I,  and  are  plotted 
in  Figure  21*  The  slopes  have  been  constructed  for  a  one  second  time  increment 
between  1.3  and  2.3  seconds.  As  long  as  this  increment  is  the  same  and  is  chosen 
at  the  same  time,  the  ratio  of  slopes,  tan  _  a  ,  remains  constant,  as 

tan  9  b 

indicated  by  Equation  (2S1.  The  time  increment  can  be  selected  between  0  and 
a  time  tj ,  of  course,  but  this  requires  knowledge  of  the  location  of  the  history 
curve's  vertex,  which  is  difficult  to  determine  experimentally  due  to  the 
response  time  of  the  system.  Selection  of  the  time  increment  is  not  entirely 
arbitrary,  however,  since  consideration  must  be  given  to  the  criterion  that 
the  solid  anpear  semi-infinite  in  thickness  during  heating,  as  satisfied  by 
Equation(20).  This  thickness  limitation,  ,  has  been  calculated  from  published 
data  at  a  time  equal  to  2  seconds  and  is  listed  in  Table  II  for  each  of  the  solids 
studied.  It  is  seen  that  specimens  having  approximately  a  one-inch  thickness, 
which  is  experimentally  convenient,  appear  semi-infinite  during  a  2  second 
observation  period.  Observation  for  periods  greater  than  this  will  result  in 
a  departure  from  the  history  curve's  parabolic  shape  after  2  seconds. 
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Figur.  21  TIME-TEMPERATURE  HISTORIES  OF  MATERIALS 


TABLE  II. 


Material  thicknesses  for  which  specimens  are  considered 
semi-infinite  in  extent. 


Material 

-  >  * 

J  at 

(t  <  2  secs  . ) 

Thickness 

L  (inches) 

=  (3  at)1/2 

Copper 

6.18 

0.98 

Aluminum 

4.98 

0.88 

Steel  (Type  1020) 

0.78 

0.35 

Lead  (Chemical  Grade) 

1.20 

0.47 

Graphite  (Type  ATJ 
with  grain) 

5.82 

0.95 

Graphite  (Type  ATJ 

against  grain) 

00 

•"O 

T 

0.82 

A  material  is  considered  semi -infinite  to  the  flow  of  heat  when  the 
condition,  L  L  „„  =  (3  dt)*^ ,  obtains.^) 

^^The  thermal  diffusivities  used  in  this  calculation  are  obtained  from 
the  published  data  listed  in  Table  I.  Although  the  experimental 
heating  time,  t,  was  3  seconds,  the  time  of  interest  in  the  calcu¬ 
lations  is  approximately  1  to  2  seconds. 


It  must  be  remarked  that  the  family  of  curves  in  Figure  21  indicate  only 
the  relative  displacements  of  each  curve.  The  group  will  appear  steeper 
the  greater  the  incident  radiation  flux  and/or  emissivity,  and  will  appear 
flatter  the  smaller  the  flux  and/or  emissivity.  If  the  flux  is  less  than  the 
critical  quantity,  FQ,  in  Equation(35  ),  during  the  heating  period  then 
observation  of  a  heating  curve  is  not  possible. 

f  .  Cooling  Method 

Another  approach  for  obtaining  thermal  conductivity  using  a  comparative 
technique,  along  similar  lines  as  described,  requires  radiometric  observation 
of  the  solid's  surface  while  it  cools  after  heating.  This  approach  is  workable, 
but  is  not  as  practical  as  the  heating  method.  Analagous  to  Equation  (21)  for 
the  temperature  history  during  heating,  is  the  following  equation,  which 
generates  the  cooling  history  at  the  solid's  surface 

T(Q,t)-T0  =  [tl/2  _  (t  _  r)l/2  ]  t  (39) 

where  r  is  the  heat  pulse  duration.  For  solution,  t  must  be  greater  than  r  . 

If  t  -  r  ,  then  Equation(Zl)  results.  Equation  (39),  then,  generates  a 
curve  which  is  the  difference  between  two  parabolas  having  the  same  latus 
rectums  but  separated  in  time  by  r  ,  even  though  reflection  exists  during 
heating.  The  chief  advantage  of  the  cooling  method  compared  to  the  heating 
method  is  that  reflected  radiation  no  longer  requires  discrimination  in  the 
experiment.  However,  sensitivity  and  signal-to-noise  ratio  must  still  be 
considered  for  observability,  and  therefore,  the  detector  and  recorder  must 
either  be  protected  against  overshooting  from  the  reflected  radiation  component 
if  monitoring  is  continuous,  or  they  must  be  activated  immediately  following 
the  termination  of  the  heat  pulse. 

One  serious  limitation  of  the  cooling  method  is  due  to  the  requirement  that 
the  solid  must  appear  semi-infinite  during  both  the  heating  and  cooling  cycler 
up  to  the  time  observation  is  complete.  Thus,  the  heating  pulse  duration  and 
the  duration  of  observation  while  the  solid  cools  together  must  be  less  than 
the  time  for  transport  of  the  pulse  through  the  body.  For  most  good  conductors, 
this  can  only  be  accomplished  by  depositing  a  sufficient  quantity  of  energy 
uniformly  over  a  suitable  region  on  the  solid's  surface,  within  a  short  interval, 
so  that  excessive  temperatures  are  not  attained,  and  then  observing  the  cooling 
curve  up  to  the  semi -infinite  time  limitation.  This  might  be  accomplished  using 
a  laser  as  a  radiation  source,  which  has  the  advantage  of  emitting  radiation  in 
other  than  the  infrared  region,  so  that  reflection  is  not  a  problem.  Analysis 
of  the  cooling  curve  to  obtain  a  tangent  ratio,  from  which  a  conductivity  ratio 
is  determined,  will  be  similar  to  that  for  the  heating  curve.  Other  than  for 
the  reflection  problem,  the  cooling  method  offers  no  advantage  over  the  heating 
method  and,  in  general,  probably  would  prove  to  be  more  critical  to  use  and 
would  require  a  more  costly  radiation  source. 


3.  Experiment 

Measurements  were  performed  on  five  planar -oriented  type  ATJ  graphite 
specimens  exhibiting  a  range  of  densities  and  four  types  of  good  thermal 
conductors,  having  a  broad  range  of  thermal  properties  and  characteristics, 
and  densities.  These  good  conductors  were:  (1)  drawn  copper;  (2)  drawn  aluminum 
(3)  type  1020  drawn  steel;  and  (4)  chemical  grade  cast  lead.  Published  values 
of  their  properties,  from  which  calculated  quantities  of  interest  are  obtained, 
are  listed  in  Tables  I  and  II.  The  purpose  of  selecting  these  materials  was 
to  determine* the  range  of  validity  of  the  technique  for  solids  having  conductivities 
near  that  of  graphite.  Copper  was  selected  principally  to  serve  as  a  standard 
for  comparison. 

a.  Apparatus 

A  schematic  diagram  of  the  experimental  arrangement  used  is  shown 
in  Figure  21.  The  radiation  source  consisted  of  a  parallel-connected  bank 
of  3  Sylvania  Type  DFA  150  watt  projection  lamps.  These  lamps  reach 
maximum  intensity  in  approximately  80  milliseconds  and  contain  built-in 
pre-focussed  reflectors,  so  that  an  object  can  be  heated  over  a  region 
whose  area  is  proportional  to  the  distance  between  it  and  the  lamp.  The 
lamps  were  arranged  to  produce  a  nearly  homogeneous  intensity  over  an 
area  having  a  3/4-inch  diameter.  Homogeneity  was  evidenced  by  scanning 
a  roughened  surface  of  a  massive  copper  heat  sink  with  an  infrared  radiometer, 
on  which  the  lamps  were  focussed,  and  observing  the  diffusely  reflected 
infrared  radiation.  A  heat  sink  was  used  so  that  self-emitted  radiation 
could  be  neglected  during  observation. 

The  infrared  detection  instrument  used  was  a  Barnes  Engineering  Co. 
Industrial  Infrared  Radiometer,  Model  R4DI.  The  infrared  wavelength 
band  passed  by  its  germanium  filter  is  1.  8  to  25  microns  (a  peak  wavelength 
of  10  microns  for  a  Planck's  black  body  radiation  distribution  corresponds 
to  a  temperature  of  approximately  300  degree  K).  This  radiometer  has 
four  desireable  features  for  this  application:  (1)  its  response  time  is  as 
fast  as  10  milliseconds;  (2)  it  can  be  easily  focussed  to  a  circular  field 
of  view  having  a  diameter  as  small  as  1/16  inch  at  a  target-detector 
distance  of  30  inches;  (3)  its  temperature  sensitivity  was  0.5  degree  C 
in  the  manner  used;  and  (4)  it  contains  an  offset  control  which  can  suppress 
reflected  radiation  so  that  small  temperature  variations  riding  on  top  of  a 
large  signal  can  be  observed  at  its  voltage  output  terminals 

A  Sanborn  Model  152-100B  two-channel,  recorder  was  used  for  all 
measurements.  Four  principle  features  of  this  unit  are:  (1)  a  writing 
response  time  of  approximately  30  milliseconds;  (2)  a  sensitivity  of  0.001 
volts/mm  deflectipn;  (3)  a  chart  speed  variable  up  to  100  mm/sec;  and 
(4)  it  contains  a  timer-marker. 
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VARIABLE  _ 

TRANSFORMER  PROJECTION 


A  photocell  was  used  to  monitor  heat  pulse  initiation  by  observing  the 
visible  radiation  reflected  from  the  specimens.  Its  signal  was  fed 
into  the  second  channel  of  the  recorder.  Although  pulse  initiation  could 
be  determined  from  the  radiometer's  signal  when  observing  the  total 
reflected  and  self-emitted  radiation  components  on  the  full  scale  of 
the  recorder,  the  photocell's  signal  was  necessary  to  determine  pulse 
initiation  when  the  reflected  component  was  suppressed  with  the  radiometer's 
offset  control. 

b.  Specimens 

1. )  Dimensions 

For  a  pulse  duration  of  two  seconds,  all  materials  studied  appear 
semi-infinite  for  thicknesses  greater  than  1  inch.  This  is  seen  from 
the  values  of  the  critical  thicknesses,  L  ,  listed  in  Table  II.  Thus, 

no 

copper,  aluminum,  steel  and  lead  were  conveniently  fabricated  into 
one-inch  cubes.  Areal  dimensions  of  1  in.  x  1  in.  were  chosen  to 
essentially  fill  the  3/4-inch  diameter  area  undergoing  heating.  Graphite 
specimens  were  conveniently  fabricated  into  blocks  that  were  one-inch 
thick  in  the  oriented  direction  with  the  grain,  and  3/4-inch  thick  in 
the  direction  against  the  grain.  The  third  dimension  was  one  inch. 

2.  )  Coatings 

All  surface  of  the  specimens  were  coated  by  spraying  with  Krylon  flat 
black  lacquer.  Coating  thickness  and  its  uniformity  were  not  controlled, 
but  it  was  estimated  that  the  average  thickness  was  no  greater  than  10"^  cm. 
This  surface  preparation  was  necessary  to:  (1)  provide  equal  and  uniform 
emissivities;  and  (2)  increase  emissivity  to  increase  incident  radiation 
absorption  and  radiation  self  emission.  It  is  recognized  that  this  lacquer 
coating  causes  resistance  to  heat  flow,  but  it  was  found  to  exhibit  the 
highest  emissivity  for  incident  radiation  of  the  several  coatings  examined. 
Chief  among  these  was  Aquadag,  which  is  a  colloidal  suspension  of  graphite 
particles.  This  substance  as  a  coating  exhibits  less  resistance  to  heat  flow 
than  does  the  lacquer,  but  its  emissivity  is  also  less.  Its  chief  disadvantage 
is  that  it  was  found  to  be  difficult  to  apply  uniformly  and  with  good  adhesion. 
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c.  Procedure 

1. )  Boundary  Conditions 

All  specimens  were  mounted  snugly  into  a  jig  fixture  which  contained 
low  conductivity  rigid  plastic  walls.  The  front  and  rear  surfaces,  of  the 
specimens,  with  respect  to  the  incident  radiation,  were  exposed.  This 
fixture  served  to  inhibit  lateral  heat  flow  loss  as  required  by  the  boundary 
conditions  of  the  technique.  Heat  loss  from  the  rtar  surface  is  actually 
desireable,  since  it  appears  to  extend  the  semi -infiniteness  of  the  specimens 
for  good  conductors.  The  three  projection  lamps  were  focussed’onto 
the  front  surface  over  a  3/4-inch  diameter  area.  For  the  conductivity 
determinations,  the  radiometer  was  focussed  into  a  1/16-inch  diameter 
area  at  the  center  of  the  heating  area.  This  diameter  ratio  of  12  was 
considered  sufficient  to  meet  the  requirement  of  plane  wave  flow 
through  the  thickness  of  the  specimen  during  the  observation  time. 

To  insure  that  the  specimens  appeared  semi -infinite  during  the 
observation  time,  the  front  surface  was  irradiated  while  the  exposed  rear 
surface  was  simultaneously  monitored  with  the  radiometer.  The  time 
delay  between  pulse  initiation  and  detection  serves  as  a  measure  of  the 
materials  diffusivity,  and  hence  critical  thickness  None  of  the  specimens 

showed  a  rear  surface  temperature  rise  greater  than  0.5  degree  C  during 
a  pulse  duration-observation  time  of  2  seconds.  Thus,  all  specimens 
were  assumed  to  be  semi -infinite. 

2. )  Emiesivity  Measurements 

Measurements  were  performed  on  all  specimens  to  determine  their 
emissivities.  This  was  accomplished  by  pulsing  the  front  surface  of  a 
coated  specimen  and  simultaneously  observing  with  the  ratiometer  both 
the  reflected  and  self-emitted  infrared  radiation  components.  Since  the 
self-emitted  component  was  small  compared  to  the  reflected  component 
the  resulting  pulse  shape  appeared  almost  square  with  slightly  canted 
skirts.  Comparison  of  the  pulse  amplitudes  of  the  black  lacquer  coated 
specimens  with  an  aluminum  painted,  massive  aluminum  heat  sink  of 
known  infrared  emissivity  ( e ^  =0,  28)  indicated  that  the  flat  black  lacquer 
coating  has  an  =  0.85,  which  compares  favorably  with  a  reported  value 
of  0.84  The  actual  recordings  for  the  aluminum  standard  and  a  coated 

steel  specimen,  from  which  the  emissivities  have  been  determined,  are 
shown  in  Figure  22.  Included  in  this  figure  are  descriptions  of  the  pulse 
shape  characteristics  and  the  geometrical  constructions  employed. 
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3* )  Conductivity  Measurements 

For  these  measurements,  the  radiometer's  offset  control  was  used 
to  suppress  almost  all  of  the  reflected  component  of  the  radiation.  The 
radiometer  sensitivity  could  then  be  increased  to  record  full  scale  the 
voltage  (temperature)  history  of  the  self-emits'ed  component.  Pulse 
duration  and  observation  were  made  for  3  seconds  rather  than  for  2  seconds 
to  observe  any  distortion  from  the  parabolic  shape  that  might  result  from 
violation  of  the  condition  of  semi-lnfiniteness.  Slopes  were  measured 
between  approximately  1.  3  and  Z.  3  seconds  after  pulse  initiation. 

Some  typical  results  of  th^se  measurements  are  shown  for  graphite 
specimen  No.  34,  in  Figure  24  •  In  Figure  (a)  an  emissivity  recording 
is  shown  in  which  the  temperature  rise  due  to  self  emission  can  be  seen. 

A  temperature  scale  is  also  included  to  illustrate  the  absolute  magnitude 
of  the  radiation  received  by  the  radiometer.  In  Figure  ?4  (b)  and  (c), 
the  actual  recordings  of  graphite  specimen  No.  34,  against  and  with  the 
grain  direction,  are  shown,  from  which  thermal  conductivities  were 
determined.  It  is  seen  that  the  surface  temperature  rise  for  each 
direction  is  only  a  few  degrees  Centigrade. 

Shown  in  Figure  2  5  is  a  thermal  conductivity  recording  of  a  large 
lead  fabrication,  whose  front  surface  was  coated  with  black  lacquer, 
and  which  has  beei.  pulsed  for  6  seconds.  Its  temperature  rise  during 
this  period  is  approximately  9  degrees  C. 

4*  Results  and  Discussion 


The  results  of  this  study  are  summarized  in  Tables  III  and  IV.  It  is  seen  thqt 
thw  Pleasured  infrared  emissivities  of  the  flat  black  lacquer  coatings  are  within 
approximately  3  percent  of  0.88.  Differences  that  occur  maybe  significant,  since 
they  are  slightly  greater  than  system  or  measurement  errors.  As  such  they  are 
attributed  to  differences  in  the  base  surface  finishes. 

a*  Loss  Factors 

Slopes  of  the  history  curves  for  the  specimens  are  listed  in  Table  III, 
together  with  their  ratios  referenced  to  the  slope  for  copper.  Comparison 
of  these  calculated  ratios  with  the  experimental  results  for  aluminum,  steel 
and  lead  indicates  that  a  reasonably  constant  loss  factor,  F},  of  about  2.  2 
exists.  It  was  not  possible  to  calculate  this  ratio  for  the  graphite  specimens 
since  their  densities  varied  and  their  conductivities  were  unknown.  This  loss 
factor  is  believed  to  arise  in  part  from  retardation  of  the  heat  flow  into  the 
specimens  caused  by  the  thermal  resistance  of  the  lacquer  :oating.  As 
indicated  in  the  Theory  Section  of  this  report,  a  consequence  of  this  retardation 
is  to  increase  the  slope  of  the  temperature  rise  to  a  degree  that  is  dependent 
on  the  heat  transfer  characteristics  of  the  material.  The  net  effect  is  to 
decrease  the  experimental  tangent  ratio  with  respect  to  its  calculated  value. 
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Fi*w»  25  A  THERMAL  CONDUCTIVITY  RECORDING  OF  A  LARGE  LEAD  FABRICATION, 
WHOSE  FRONT  SURFACE  ISCOATEO  WITH  FLAT  BLACK  LACQUER 
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3.  3  was  assumed  for  all  graphite  apecimeni  in  the  calculations 

neasured  correspond  d  to  the  1.8  to  25  mlc-on  wavelength  bond  observed  by  the 

meter's  germanium-fi  bred.,  bolometer  detector. 
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TABLE  IV.  Densities  and  experimental  ratios  and  values  of  k  for 
Type  ATJ  graphite  samples  for  heat  flow  with  and 
against  the  direction  of  the  grain  orientation. 


Graphite 

Sample 

Number* 

Density 

P 

(gm/cu  cm) 

** 

kw 

kAg 

It//  cal  i 

lr  +  f  \ 

Ag  cm-secPC) 

w'crmsec«C) 

13 

1.700 

1.02 

0.13 

0.13 

16 

1.781 

1.11 

0.17 

0.19 

29 

1.649 

1.19 

0.13 

0.16 

31 

1.819 

1.06 

0.23 

0.24 

34 

1.600 

1.19 

0.11 

0.13 

$ 

The  dimensions  of  ail  graphite  samples  were  3/4"  x  1"  x  I".  The 

3/4"  x  1"  face  was  perpendicular  to  the  heat  flow  direction  with  the  grain. 

For  highly  oriented  type  ATJ  graphite  the  quantity  kw/k^g  is  calculated 
to  equal  1.33  from  published  data. 

+A  loss  factor  of  3.3  was  used  to  obtain  these  k's  for  graphite. 


Another  contributing  effect  to  F1  is  attributed  to  later nal  heat  flow 
loss  which  results  in  nonplanar  wave  pulse  transmission  through  the 
specimen.  Lateral  flow  makes  a  specimen  appear  to  be  a  better  conductor 
than  it  really  is.  This  loss  is  more  pronounced  in  a  poorer  conductor 
so  that  its  net  effect  is  to  reduce  the  experimental  tangent  ratio  also,  thus 
increasing  the  loss  factor,  Fj,  further.  The  fact,  however,  that  Fj  appears 
nearly  constant,  which  indicates  a  constant  shift  in  experimental  values 
from  their  calculated  values  does  provide  a  means  for  easily  adjusting 
this  difference 

The  {k  p  Cp  )  ratios  referenced  to  copper,  from  which  the  conductivities 
are  determined,  are  also  included  in  Table  III,  together  with  their 
corresponding  calculated-to-experimental  loss  factors,  F^.  These 
factors  are  seen  to  center  about  a  value  of  3.  3,  and  are  similar  in  their 
origins  as  for  Fj. 

b.  Graphite  Conductivities 

Listed  in  Table  IV  are  the  densities  and  the  experimentally  determined, 
grain-oriented  conductivity  ratios  for  the  5  graphite  specimens  studied. 

It  is  seen  that  the  conductivities  in  a  direction  with  the  grain  are  consistently 
greater  than  the  conductivities  against  the  grain,  regardless  of  density. 

These  ratios  are  to  be  compared  to  a  value  of  1.33  calculated  from  publishe^ 
data  for  highly  oriented,  Type  ATJ  graphite  having  a  density  of  1.  7  gms/cm  . 

In  partial  explanation  of  this  difference,  it  is  to  be  noted  that  the  experimental 
technique  necessarily  results  in  values,  which  are  representative  of  only  a 
localized  region,  whereas  the  densities  of  the  specimens  are  bulk  values. 

Also  included  in  Table  IV  are  the  conductivities  of  the  graphite  specimens 
in  directions  with  and  against  the  grain.  These  values  were  determined  using 
a  loss  factor,  F^,  equal  to  3.3  and  referenced  to  copper.  These  experimental 
conductivities  are  to  be  compared  to  published  values  of  0.  28  and  0.  21 
cal/cm-sec-degC,  with  and  against  the  grain,  respectively.  The  experimental 
values  are  seen  to  be  generally  lower  than  the  reported  values,  which  suggests 
that  either  F2  is  less  than  3.  3  for  graphites,  or  the  specimens  have  lower 
conductivities  in  actuality  than  the  reported  values.  Unfortunately,  steady-state 
conductivity  measurements  were  not  performed  on  any  of  the  specimens  due 
to  their  geometry.  However,  had  these  measurements  been  performed,  they 
would  have  validity  only  for  the  bulk  specimens.  The  experimental  conductivities 
of  the  specimens  as  a  function  of  their  densities  are  plotted  in  Figure  26. 

It  is  seen  that  conductivity  is  proportional  to  density,  as  expected,  such  that 
a  one  percent  change  in  density  produces  approximately  a  6  percent  change 
in  conductivity. 
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c.  Error  Analysis 

Contributions  to  the  errors  involved  in  determining  the  thermal  conductivities 
arise  from:  (1)  instrumentation;  (2)  geometrical  construction  on  the  history  curves; 
and  (3)  differences  between  published,  or  "standards",  values  and  experimental 
values.  Errors  in  the  precision  of  the  instrumentation  are  estimated  to  be  within 
2  percent.  Since  a  comparative  technique  was  used,  knowledge  of  the  errors  in  the 
accuracy  of  the  instrumentation  is  unnecessary. 

The  probable  error  in  determining  thermal  conductivity  is  estimated  to  be 
17  percent.  This  estimate  was  arrived  at  by  determining  the  error  from  each 
contributing  factor  in  Equation  (32)  and  then  calculating  the  probable  error  by 
taking  0.67  of  the  square  root  of  the  sum  of  the  squares  of  each  error  contribution. 

The  emissivity  error  was  estimated  as  2  percent  and  the  maximum  slope  error  as 
21  percent.  The  slope  error  arises  from  the  uncertainty  in  constructing  the  mean 
of  a  history  curve  when  system  noise  is  significant.  This  uncertainty  can  be 
diminished  greatly  if:(a)  a  longer  observation  time  is  used  to  increase  the  abscissa 
of  the  slope  (which  requires  using  larger  specimens);  and/or  (b)  a  more  intense  heat 
source  is  used  to  increase  the  slope  of  the  temperature  rise  and  suppress  noise  by 
reducing  radiometer  sensitivity.  If,  for  example,  a  construction  time  of  2  seconds 
had  been  used  instead  of  the  1  second  selected  in  this  study,  the  probable  e.ror  in 
determining  conductivity  would  have  been  approximately  9  percent  rather  than  17  percent. 

The  scatter  in  the  conductivity  vs.  density  data,  plotted  in  Figure  26,  is  attributed 
mainly  to  the  construction  error.  An  additional  source  of  contribution  to  the  scatter 
is  attributed  to  local  density  variations  existing  within  a  specimen,  having  a  measured 
bulk  density,  which  influences  the  history  curve  during  the  observation  period. 
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V. 


CONCLUSIONS  AND  RECOMMENDATIONS 


Highlights  of  the  capability  developed  during  this  contract  are  summarized 
below.  It  is  to  be  noted  that  many  of  these  conclusions  were  identified  during 
the  first  year's  work,  and  have  been  extended  to  other  grades  of  interest  to 
the  Air  Force  during  the  second  year. 

1)  The  predominant  material  variable  influencing  the  properties  of 
interest  in  bulk  graphite  for  ablative  applications  has  been  identified 
and  experimentally  verified  to  be  bulk  density.  Total  volume  porosity, 
which  can  be  calculated  from  bulk  density  measurements  via  graphite's 
theoretical  density,  of  discrete  volumes  has  been  noted  as  an  important 
parameter  in  studies  of  erosion  rate  in  rocket  nozzle  inserts. 

2)  A  state-of-the-art  nondestructive  radiation  gauging  technique, 
which  was  developed  during  the  first  year,  has  been  extended  to 
all  grades  of  interest,  and  can  be  used  to  determine  density  in 
discrete  volumes  cf  billets  and  hardware  of  these  grades  to  within 
±1%.  As  mentioned  above,  this  capability  leads  to  the  determination 
of  total  volume  porosity  within  the  same  discrete  ’’■olumes  where 
radiation  measurements  are  taken. 

3)  Nondestructive  ultrasonic  longitudinal -wave  velocity  measure¬ 
ments  have  been  correlated  to  elastic  modulus,  both  with  and 
against  the  grain,  for  the  grades  investigated.  Combined  use  of 
velocity  measurements  and  radiometrically  determined  density 
have  been  correlated  uniquely  to  elastic  modulus  to  provide 
determination  of  this  property  to  withiniO.06  x  10^  psi  and,  by 
virtue  of  a  tensile  test  parameter  identified  herein  as  the  secant 
modulus,  such  measurements  can  also  be  used  to  determine 
ultimate  tensile  strength  to  within  ±300  psi.  It  is  particularly 
significant  to  note  that  tripling  the  range  of  values  of  these 
mechanical  properties  correlations  to  include  new  grades  studied 
within  the  past  year  does  not  alter  the  correlations  significantly. 

4)  Based  upon  the  thermoelastic  experimentation  and  analysis  made 
during  this  contrast,  as  well  as  the  related  studies  of  wave  propagation 
in  anisotropic  elastic  media,  it  seems  plausible  that  a  qualitative 
NDT  parameter  for  comparing  the  relative  resistance  of  graphite 
materials  to  thermal  shock  can  be  obtained  empirically.  Verification, 
of  course,  will  require  extended  study. 
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5)  Based  on  limited  experimentation  reported  herein,  the  determination 
of  the  thermal  conductivity  of  small  graphite  fabrications  near  room 
temperature  is  feasible  using  a  nondestructive,  one  surface,  transient 
heating,  comparative  technique.  It  is  estimated  that  the  accuracy  of  this 
infrared  technique  is  pnesaitly  within  177o  of  the  standard  used  for  comparison. 
Extension  of  this  technique  to  include  large  fabrications  is  also  feasible. 
Furthermore,  it  is  believed  that  the  technique  can  be  applied  to  solids, 

in  general,  at  elevated  temperatures.  Validation  of  the  technique,  however 
will  require  extended  study. 

6)  Using  this  technique  it  was  noted  that  a  one  percent  change  in  the 
density  of  ATJ  graphite  corresponds  to  as  much  as  a  6  percent 
change  in  its  thermal  conductivity,  illustrating  the  influence  of 
density  on  a  thermal  property,  and  the  sensitivity  of  this  technique 
to  variations  of  this  magnitude.  Limitations  in  the  accuracy  for 
nondestructive  infrared  measurement  of  conductivity  are  those 
associated  with  the  mechanics  of  geometrical  construction  used  in 
the  technique. 

7)  Practical  aspects  of  the  accomplishments  noted  above  are  that, 
within  the  limits  of  the  NDT/properties  correlations  described, 
determination  of  density/porosity  and  mechanical  properties  can 
be  made  by  use  of  these  state-of-the-art  techniques  in  discrete 
volumes  on  billets  and,  with  minor  technique  adjustments  for 
changes  in  geometry,  on  finished  hardware  components,  thereby 
providing  an  assessment  of  these  properties  for  increased  reliability, 
which  would  otherwise  not  be  possible.  The  feasibility  of  a  practical, 
one-sided,  infrared  technique  for  measuring  thermal  conductivity 
has  also  been  demonstrated. 

The  following  recommendations  are  offered  on  the  basis  of  progress 
made  to  date  and  described  in  this  report,  and  the  knowledge  obtained 
through  these  studies  aj  to  specific  areas  where  consideration  should 
be  given  to  extended  and  additional  studies. 

1)  It  is  recommended  that  the  correlation  program  and  NDT  technique 
development  be  continued  on  graphite  covering  the  individual  and 
combined  application  of  radiometric,  ultrasonic,  and  infrared 
tests  to  obtain  correlations  of  a  unique  nature  with  properties  and 
behavior  characteristics  of  interest  for  ablative  applications. 

i)  It  is  recommended  that  better  insight  into  thermal  shock 
phenomena,  as  well  as  failure  mechanisms  in  general,  can 
be  achieved  by: 

a.  experimental  investigation  of  wave  phenomena  at 

elevated  temperatures 

b.  triaxial  tests  over  a  range  of  temperatures  and 

o  nentations 

c.  routine  tension,  compression  and  torsion  tests  at 

orientation  intermediate  to  "against"  and  "with 
the  grain"  directions 

d.  multiaxial  stress  experiments 
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3)  It  is  recommended  that  infrared  technique  development  and  related 
studies  be  extended  to  determine  limitations  that  may  exist,  and  to 
improve  on  these  where  possible,  as  well  as  to  adapt  these  techniques 
to  practical  field  use.  Specifically,  work  is  necessary  to: 

a.  determine  the  effect  of  coatings  and  their  thickness  on 

retardation  of  heat  flow  into  the  solid 

b.  establish  the  uniformity  of  heat  flux  and  the  ratio  of 

diameters  of  the  heating  area  and  radiometer  field  of 
view  necessary  to  achieve  a  desired  degree  of 
reliability  and  accuracy  in  fabrications  having 
various  configurations 

c.  determine  the  radiation  source  strength  and  type  most 

suitable  for  various  applications 

d.  determine  minimum  instrumentation  requirements 

e.  explore  methods  to  facilitate  data  readout  and  analysis  <■ 

4)  It  is  recommended  that  action  be  taken  to  promote,  demonstrate 
and  field  test  the  applicable  NDT  for  benefit  of  producers,  fabricators, 
and  users  of  graphite  materials.  Such  activity  will  assist  in  the 
development  and  qualification  inspection  of  improved  graphite  grades, 
and  better  identify  the  failure  modes  and  mechanisms  which  frustrate 
the  development  of  design  models  for  ablative  applications. 
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Appendix  I. 


Wave  Propagation  in  Transversely  Isotropic  Media 
a.  Wave  Propagation 

The  stress -strain  laws  of  plane  isotropy  can  be  written  as 
°  1  v2\ 


<1 

<2 

*3 

Y\2 


1 


1 


(a2  +  °3  ) 


'2\ 


v2\ 


v2  a3 


c2 

v2  a2 


12 
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These  stress-strain  laws  can  be  used  with  the  strain  displacement  and  the  equilibn  :m 
equations: 
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Kolsky  ,  as  well  as  many  others,  has  presented  the  correct  wave  equations 
for  anisotropic  elastic  media.  For  the  sake  of  completeness  of  discussion, 
they  are  listed  here.  As  an  example  of  solutions  of  the  equations  of  motion, 
consider  the  case  of  plane  waves  in  an  infinite  transversely  isotropic  medium. 
Take  a  wave  with  its  normal  in  the  direction  of  the  axis,  X^  . 


U;  *  Aj  cos  (  X,  exp 


j  =  1,  2,  3 


These  may  now  be  substituted  into  the  displacement  equation  of  motion. 

The  corresponding  secular  equation  can  then  be  derived.  After  applying 
the  necessary  condition  for  non-trival  solution,  i.e. ,  setting  the  determinant 
of  the  equation  equal  to  zero: 
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with  the  resulting  wave  velocities 
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For  graphite,  because  of  the  small  Poisson's  ratio  2s  .06,  at  .10) 
for  room  temperature,  the  wave  velocity  can  be  very  closely  approximated 

by  '/Ej/p. 

V 

However,  for  elevated  temperatures,  these  Poission's  ratios  increase  to  .26 
and  .40.  In  addition  to  the  marked  thermal  effects,  apparent  directionally 
dependent  strain  rate  influences  have  been  reported^',  for  the  Poisson's  ratio. 

b.  Coupled  Thermoelastic  Stress  Wave  Phenomena 

Existing  thermal  stress  analysis  for  anisotropic  elastic  media  is  based 
on  linear,  uncoupled,  quasi-static  thermoelastically  and  therefore  the 
solutions  are  satisfactory  for  equilibrium  conditions^®'  ^9).  Weng, 
for  example,  has  considered  the  problem  of  thermal  stresses  in  anisotropic 
hollow  cylinders  and  has  shown  some  of  the  effects  of  anisotropy.  He  reports 
thermal  shock  parameters  of  the  form 
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However,  in  the  strictest  sense,  no  reasonable  mathematical  solutions 
have  been  obtained  for  the  thermal  shock  problem  involving  the  determination 
of  the  transient  temperature  and  stress  distribution. 

By  appealing  to  the  appropriate  stress -strain  laws,  the  equilibrium 
equations,  conservation  of  mass  and  energy  and  the  laws  of  thermodynamics, 
the  coupled  anisotropic  thermoelastic  field  equations  can  be  formulated, 

In  particular,  the  partial  differential  equation  governing  the  temper ature 
field  is 
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where 
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K-  i  =1,2  are  thermal  conductivities 

*  l 

■i  ,  ^  =  thermal  expansion  coefficients 

The  equation  is  unmanageable  but  serves  to  illustrate  the  complexities 
introduced  by  the  more  realistic  field  equations. 


It  is  interesting  to  compare  equation  (16)  to  the  corrected  "heat"  equation 
of  isotropic  thermoelastic  medial®*  ^1). 


The  directionally  dependent  materials  properties  of  (16)  are  seen  to 
modify  the  transient  thermal  response.  Engineering  calculation  of 
thermal  stresses  are  ordinarily  carried  out  in  two  steps: 


1)  an  estimate  of  the  temperature  distribution  is  obtained  from 
Fourier's  heat  conduction  equation 


=  ?CvT 


2)  next,  the  stresses  are  calculated  from  the  equations  of 
elasto-statics. 


Comparing  this  procedure  to  the  solution  of  the  more  exact  equation  (17), 
one  can  see  that  two  effects  are  omitted:  the  effect  on  temperature  distribution 
of  straining  the  body  and  the  dynamic  effect  due  to  inertia  forces. 
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Appendix  11:  Machanical  Propartiaa 
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♦Object  in  selecting  specimens  from  each  billet  was  to  obtain  the  maximum  range  in  these  properties  as  indicated 
by  the  NDT  response  rather  than  to  obtain  a  statistical  average  of  random  samples. 


APPENDIX  IQ. 


Test  Methods  for  Determination  of  Physical  Properties  of  Graphite 


1.  Bulk  Density:  Bulk  density  was  determined  by  dividing 
the  mass  of  a  specimen  by  its  bulk  volume.  The  conventional 
unit  is  grams /cc 

2.  Specific  Resistance:  Specific  resistance  was  determined 
by  measuring  the  potential  drop  along  the  length  of  a  specimen 
of  known  geometry  with  a  given  current  flow.  Measurements 
were  made  at  room  temperature. 

3.  Dynamic  Modulus:  Room  temperature  modulus  measure¬ 
ments  were  performed  by  the  resonant  frequency  vibration 
technique  with  the  specimen  supported  at  two  points  as  a  beam. 

a.  Equipment:  Type  FM-500  Elastomat,  Magna  flux  Corp. 

b.  Specimen  Geometry:  1 1/8"  diameter  by  6"  nominal  length 

4.  Tensile  Properties 

a.  Equipment:  Instron  Model  TTC  tensile  testing  machine 

b.  Strain  Rate:  Constant  to  failure  0.005  in/in/min. 

c.  Temperature:  75*  F  ±2*F 

d.  Tensile  Specimen:  Gver-all  length  of  5  inches,  gauge 
diameter  0.400  inches. 

e.  Strain:  Measured  by  ASTM  Class  B  1  extensometers 
having  a  2-inch  gauge  length 
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